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PHYSICS  AI^D 
MTHEMATICS 

PREFACE 

We  use  the  term  pre-Maxwell  equations  to  signify  a 
formalization  of  the  electromagnetic  theory  which  existed 
prior  to  Maxwell's  introduction  of  the  concept  of 
displacement  current.   The  purpose  of  this  report  is  to 
give  a  general  descriptive  survey  of  the  theory  of  these 
equations  which  is  aimed  primarily  at  applications.   The 
pre-Maxwell  equations  yield  satisfactory  approximations 
to  solutions  of  electromagnetic  problems  in  a  niomber  of 
fields;  chiefly,  in  classical  network  and  electrical 
machinery  problems,  and  in  large  parts  of  fluid  magnetics. 
The  elimination  of  electromagnetic  wave  propagation  at 
the  very  outset  in  problems  where  it  does  not  play  an 
Important  role  frequently  results  in  great  simplification. 
Moreover,  in  formulations  of  fluid  magnetics  which  use 
non-relativistic  fluid  dynamic  equations,  the  pre-Maxwell 
equations,  being  Galilean  invariant,  are  in  some  respects 
even  more  appropriate  than  Maxwell' s  equations. 

The  system  of  pre-Maxwell  equations  has  its  own 
intrinsic  mathematical  interest  because  of  properties 
which  are  entirely  different  from  those  of  Maxwell's 
equations.   It  is  possible  to  give  a  rather  complete 
existence  theory  for  the  system  of  pre-Maxwell  equations 
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by  using  the  properties  of  harmonic  vectors.  Many  of 

the  results  quoted  here  have  been  proved  in  recent 

+ 
years  ,  however  we  have  not  hesitated  to  label  as 

"theorem"  a  number  of  additional  results  which,  although 

not  yet  rigorously  proved,  seem  susceptible  of  proof 

by  the  same  techniques. 

The  valiiable  assistance  of  Hanan  Rubin  in  preparing 

an  early  draft  of  this  report  is  gratefully  acknowledged. 

The  fundamentals  of  this  theory  were  originally  outlined 

in  195U. 


H.  G. 


•f  For  references  see,  e.g.,  K.  0.  Priedrichs,  "Differ- 
ential forms  on  Riemannian  manifolds,"  Comm.-on  Pure 
and  Appl.  Math.,  VIII,  551-590  (1955). 
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NOTES  ON  M&ONET0-HY3>R0DYNAMICS  -  NUMBER  V 
Theory  of  Maxwell's  Equations  without 
Displacement  Current 

A.  A.  Blank,  K.  0.  Priedrichs,  H.  Grad 

1»   Introjduction 
In  many  applications  of  the  electromagnetic  equations, 
e»g,  to  fluid  magnetics,  the  displacement  current 
term  in  Majcwell'  s  equations  is  considered  negligibly 
small  and  therefore  is  omitted.   The  truncated  electrQ- 
magnetic  system  , 

i|  +  curl  E  =  0  div  E  =  ^ 

(1.1) 

curl  B  =  p-  J  div  B  =  0 

obtained  in  this  way  will  be  called  the  pre-Maxwell  system 
because  it  represents  the  content  of  electromagnetic  theory 
before  the  introduction  of  displacement  current  by  Maxwell. 

These  equations  are  an  appropriate  formulation  of 
electromagnetics  in  the  conventional  engineering  analysis 
of  rotating  machinery.   In  fluid  magnetics,  when  the  fluid 
equations  are  given  in  non-relativistic  foi»m,  the  pre- 

Maxwell  system  would  seem  the  more  appropriate  for  its 

2 
invariance  under  Galilean  transfoiroation  ,  compared  with 

1.   In  rationalized  units* 

2«   Galilean  invariance  is  obtained  by  transforming  as 

vectors  the  quantities  Eu  =  E  +  u  >*  B   (polar)  and  B 
(axial).   Eu   and  B   are  to  be  interpreted  as  the 
electric  and  magnetic  field  intensities  as  observed  in 
the  given  frame  moving  with  velocity  u» 
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the  Lorentz  invarlance  of  Maxwell's  equations.  This 
usage  will  not  only  prove  less  awkward  but  may  possibly 
be  more  accurate. 

It  is  usually  assumed  that  the  equations  (1.1) 
represent  an  approximate  description  of  a  certain  class 
of  solutions  of  the  complete  system  of  Maxwell's  equations. 
Actually,  although  the  difference  in  the  two  systems  of 
equations  may  appear  slight,  it  cannot  be  assumed  without 
analytical  investigation  that  a  solution  of  one  will 
approximate  a  solution  of  the  other.   In  fact,  despite 
the  similarity  of  the  two  systems,  and  even  though  it  is 
true  that  the  solutions  of  one  may  under  certain  circvim- 
stances  approximate  solutions  of  the  other,  the  mathematioal 
properties  of  the  two  systems  differ  radically.  The  differ- 
ence between  the  two  systems  is  closely  comparable  to  the 
difference  between  Laplace's  equation  and  the  wave  equation. 
For  example,  the  pre-Maxwell  system  admits  essentially 
only  boundary  values,  where  Maxwell's  system  admits  not 
only  boiindary  values,  but  initial  values  also. 

To  gloss  over  the  differences  between  the  two  systems 
may  be  serious.   For  example,  the  system  (1.1)  is  sometimes 
used  together  with  the  equation  of  conservation  of  chaj?ge 
(^q/()t  =  div  J,   despite  the  fact  that  this  relation  will 
usually  contradict  the  equation   div  J  =  0,   a  plain 
consequence  of  the  equations  (1,1  )<,   The  pre-Maxwell  system 
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together  with  charge  conservation  is  overdetennined  and 
will  in  general  have  no  solution  at  all.  This,  of  course, 
does  not  preclude  the  possibility  that  a  given  solution 
of  the  pre-Maxwell  system  may  be  a  very  good  approximation 
to  some  solution  of  Maxwell's  equations  which  does  satisfy 
conservation  of  charge.   For  a  solution  of  (1,1)  which  is 
a  satisfactory  approximation  to  a  solution  of  Maxwell's 
equations,  it  is  to  be  expected  that  a  slight  additional 
current  would  suffice  to  give  conservation  of  charge. 
Nonetheless,  the  important  mathematical  point  is  that 
conservation  of  charge  must  not  be  imposed  upon  the  system 
(1,1),  for  it  will  generally  create  an  insoluble  problem, 

A  crucial  question  which  should  be  asked  in  connection 
with  a  given  system  of  equations  is  whether  a  solution 
exists  corresponding  to  every  physically  relevent  problem. 
Equally  Important  is  the  knowledge  that  there  is  no 
ambiguity,  that  the  conditions  of  a  problem  determine  its 
solution  uniquely.   If  a  stated  problem  admits  precisely 
one  solution  it  is  said  to  be  well-posed,  A  well-posed 
problem  for  the  Maxwell  system  is  of  a  character  rather 
different  from  that  of  a  well-posed  problem  for  the  pre- 
Maxwell  system.   For  the  system  (1,1)  a  well-posed  problem 
is,  except  for  slight  modifications,  essentially  a  pure 
boundary  value  problem  despite  the  appearanoo  of  the  time 
derivative.   We  shall  see  that  initial  values  in  a  pre- 
Maxwell  problem  may  be  assigned  to  certain  discrete 
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parameters   and  to     E     or     B     upon  the  boundary  only,   but 
not   to     E     or     B     over  the  entire   domalno      The  mathemati- 
cal distinction  between  the    systems   Is   that   the   Maxwell 
system  Is   of  hyperbolic    type  while   the  pre-Maxwell   system 
Is  almost  elliptic   In  character.      The  difference  Is 
exemplified  by  the  fact   that  no  property  akin  to  wave 
propagation  may  be   attributed  to    the   solutions   of   (l«l)« 

As   the    simplest   Illustration  of  these  differences, 
we  present  the  following  typical  well-posed  problems  for 
the   two   sets  of  equations*     In  each  case  we  take     q  «=  0, 
J  =  0. 

a)e     The  domain  Is   the  Interior  of  a   sphere.     For 
both  systems   of  equations,   the  normal   components,   En.     and 
Bn,      are   specified  arbitrarily     as  functions  of  time  on 
the    surface     S     of   the   sphere      subject  only   to   the   condi- 
tions   /Ends  =    /B^dS  =  Oe      In   the   Maxwell   system,   initial 
data     E     and     B     are   given  at      t  =  0     subject   to   the 
conditions      dlv  E  =  div  B  =  0,      In  the   pre-Maxwell   system 
no    initial   data  can  be    imposed  at    allf   the   solution  of 
(lol)   is   already  uniquely  determined  as   a  function  of   time 
by  the  boundary  conditions « 

b).      The  domain  is   again  the   interior  of  a   sphere. 
In  this   case,   however,    the   tangential   component     E-j;     of     E, 
is   specified  on  the    surface   of  the   sphere    as    a  fionction  of 
time.      For  Maxwell's  equations,   we  would  further  specify 
E     and     B      initially  subject   to      div  E  =  dlv  B  =  0, 
However,    to  make   the   pre-Maxwell  problem  well-posed,   we 
may  further    specify  only  the   initial   values  of     Bq     on 
the   s^lrface   of   the   sphere    (subject   to    /B^dS  =  0)j 
nothing  can  be   initially  prescribed  over   the  whole 
domain© 

c).      Let   the  domain  be   the   interior  of  a  torus. 
For  the  Maxwell   equations   the    situation  is   exactly  the 
same  as   in   the   preceding  example   for   the    sphere |   we   can 
prescribe   the  boundary  condition     Et  =  0     as   well   as 
initial   values  of     E     and     B      Inside,    restricted  by 
div  E  =  dlv  B   =  0.      For   the  pre-Maxwell  equations,    in 
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addition  to  the  boundary  condition  Et  =  0  for  all  time, 

we  can  specify  B^^  =  0  initially,  just  as  in  the  previous 

example,  and  also  the  initial  value  of  the  flux  JB  •  dS 
computed  on  a  disc  which  cuts  the  torus« 

While  the  Maxwell  equations  admit  essentially 
arbitrary  initial  conditions,  we  see  that  the  pre-Maxwell 
equations  yield  a  unique  solution  tinder  much  less  stringent 
initial  conditions}  in  the  above  examples,  these  conditions 
involved  the  initial  prescription  of  boundary  values,  or 
the  value  of  a  single  parameter  (the  magnetic  flux),  and 
in  one  case  there  was  no  initial  condition  r.t  all.   The 
general  situation  is  described  in  the  theorems  of  Section  10, 

The  two  sets  of  equations  are  further  ^.istinguished  by 
the  fact  that  the  appropriate  singular  solutions  of  (1«1) 
have  a  character  distinctly  different  from  those  of  the 
Maxwell  system.   Corresponding  to  an  electron  at  rest, 
there  is  a  solution  of  either  system, 

E  =  grad  (^)   ,    B  =  0   . 

For  Maxwell's  equations  there  also  exists  a  singular 
solution  corresponding  to  a  moving  point  charge*   There 
is  no  corresponding  solution  of  (1,1)  in  which  the  moving 
charge  is  considered  to  be  a  localized  singular  current 
source  as  well  as  a  charge*  This  is  clear,  since  (lol) 
implies   div  J  =  0.   In  other  words,  the  physical  pictxxre 
of  discrete  moving  charges  must  be  smeared  into  an 
incompressible  current  flow  before  approaching  the  system 
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(1,1).      A  symptom  of   this   state  of  affairs   is    the  fact 
that   the  usual   formula  given  for  the  magnetic  force 
between  two   current   elements    (with  no  retardation)   does 
not  satisfy  Newton's   third  law   J  note   that   this  is  not 
a  consequence  of  the  neglect  of  radiation  moraent\mi»      If 
the   current    elements  are   completed  into   closed  loops 
(corresponding  to  the   singular  form  of  the   condition 
div  J  =  0),      the  theory  of   (1.1),   and  consequently  of 
this  force  formula,    is  applicable  and  yields  equal  action 
and  reaction. 

Tizming  again  to   the  general   theory  of   the  pre- 
Maxwell   system   (1.1),  we  note  that   the  form  of  the   system 
suggests   a  certain  simplification J   namely,    that   a  problem 
may  be   solved   in  two  parts,   first   for     B     and  then  for     E, 
or  vice  versa.      For   example,    we  may  begin  by  obtaining     B 
first  from   the    system 

1  CTorl   B  =  LL J 

(1.2)  < 

1     div   B   =  0 

Then,  utilizing  the  values  of  B  obtained  from  (1.2),  we 
may  determine  E  from 

rcurl  E  =  -  ^ 

(1.3)  < 

(^  div  B  =  1^ 

^     E.g.,    StT«atton,    J.   A.,   Electromagnetic  Theory, 
McGraw-Hill,,   I9I4.I,   Pa  26fe,   Ex.    22, 
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The  systems  (1.2)  and  (1«3)  are  closely  related  to  the 
potential  equation.   In  fact,  if  J  =  0,   B  has  a  scalar 
potential  which  satisfies  the  Laplace  equation.  The 
inhomogeneous  equations  (1«2)  and  (1.3)  have  a  theory 
very  similar  to  that  of  Laplace's  equation.   In  particular, 
the  natural  problems  are  boimdary  value  problems*  In 
(1»2)  if  J  is  assumed  given  as  a  function  of  position 
and  time,   B  may  be  fixed  uniquely  by  boimdary  data  and 
values  of  certain  parameters  called  periods  which  later 
will  receive  their  physical  Interpretation  as  currents 
or  fliixes^»   Once  B  is  known,  the  vector  E  may  then 
be  found  subject  to  the  prescription  of  similar  boundary 
conditions  and  periods.   In  this  case  the  periods  are 
naturally  interpreted  as  electric  charge  or  as  e,ra.f»  • 

Alternatively,  it  Is  possible  to  solve  (1.1)  by 
first  solving  for  E  In  the  system, 

curl  curl  E  =  -jx  4r 
div  E  =  I 

obtained  by  eliminating  ^B/^t     from  the  system  (1«1). 
If  E  is  determined  uniquely  as  the  solution  of  a  well- 
posed  problem  for  this  system,  then,  without  the 

4»  A  ciiri^ent  is  defined  by  the  period  X/^fQ  •  dx  defined 
on  a  closed  curve |  a  flux  by  /b  •  dS  over  a  diac» 
These  concepts  will  later  be  made  more  precise* 

5«   The  total  charge  within  a  surface  is  given  by 

ic^E  •  dS.   The  e.m.f.  or  electromotive  force  along 
an  arc  is  given  by  /E»  dx. 
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prescription  of   additional  boxmdary  data,      B     may  be 
foTjnd  from 

(1.5)  ^  "   "  ""^^^  ^      * 

only  certain  Initial  properties  of  B  will  be  needed* 
In  this  problem  all  the  boundary  data  are  Imposed  upon 
the  vector  B  In  contrast  with  the  preceding  problem 
where  boundary  conditions  are  divided  between  the  two. 

A  well-posed  problem  for  the  system  (1»1)  will 
generally  Involve  boundary  and  period  data  and  initial 
values  of  such  data.  "We  have  seen  that  the  data 
appropriate  to  the  method  of  solving  sequentially  {1.2) 
and  (1.3)  are  not  exactly  the  same  as  the  data  appropri- 
ate to  the  sequence  (Li^.),  (1.5).   Roughly,  the  difference 
lies  in  the  way  the  subsidiary  data  have  been  apportioned 
to  E  and  B.   To  formulate  the  general  well-posed  problem 
for  (1.1)  without  reference  to  a  particular  method  of 
solution,  we  shall  characterize  compatible  sets  of 
boiindary,  period,  and  initial  data  for  both  E  and  B, 

In  more  complicated  cases,  as  in  fluid  magnetics, 
J  and  q   are  not  given  a  priori  but  are  to  be  deter- 
mined together  with  B  and  B   after  supplementing  (1.1) 
by  further  equations.   For  example.  Ohm's  law  connects   J 
with  E   and  B   as  well  as  the  motion  of  the  medium* 
Compatibility  conditions  between  B   and  E  may,  of 
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coiirse,  be  modified  in  such  cases.   Frequently,  it  will 
turn  out  that  the  relation  dlv  E  =  q/K   is  not  used  in 
determining  E,  but  only  to  obtain  q  after  E  has 
been  found  otherwise.  This  situation  will  be  found  to 
obtain  for  a  perfectly  conducting  fluid,  a  case  which 
will  be  discussed  at  length  in  subsequent  notes. 

The  study  of  the  pre-Maxwell  system  is  reduced  by 
(1.2)  and  (1.3)  to  the  study  of  simpler  systems  having 
the  form 

r  dlv  Z  =  «r 

(1.6)     1         _  ,    _    , 

l^curl  Z  =  j:i  Idiv  j^  =  0). 

As  is  often  the  case  In  linear  analysis,  many  '^f  the 
difficulties  in  the  solution  of  the  inhomogeneous  system 

(1.6)  may  be  assigned  to  the  solution  of  the  homogeneous 
system 

f  div  Z  =  0 

(1.7)  i 

(^  curl  Z  =  0   . 

The    solutions   of   system   (1.7)   are   the    so-called  harmonic 
vectors. 

Harmonic   fields  may  arise  natxorally   In  several  ways. 
If      J  =  0     in   some   domain  then     B      is  harmonic   there.      If 
q  =  0      and     B      is   independent   of   time   then     E     is  harmonic. 
If  the  domain   is   simply  connected,    then  a  harmonic     B     is 
the   gradient   of   a  scalar  potential   and  the  problem  of 
determining     B     reduces   to   that  of   solving  the   Laplace 

-  18  - 


equation.   This  approach  to  the  problem  is  quite  general 

and  any  specific  additional  complications  can  only  arise 

from  the  possible  multi-valuedness  of  the  potentials, 

essentially  a  topological  difficulty.   In  applications 

to  fluid  magnetics  where  we  are  concerned  with  the 

magnetic  field  in  a  vacuum  domain  bounded  by  conductors, 

it  is  common  to  impose  the  normal  boundary  condition 

B  =0.  Vlhen  B  is  harmonic  this  condition  implies 
n 

that  B  must  vanish  identically  in  any  simply-connected 
domain  «   However,  if  the  vacuum  domain  is  bounded  by  a 
torus,  this  conclusion  does  not  follow;  there  does  exist 
a  regular  non-trivial  field.   For  a  precise  analysis  of 
these  problems  it  is  necessary  to  introduce  a  few  topologi- 
cal concepts  vftiich  clarify  the  pertinent  intuitive  notions. 

In  the  first  few  sections,  then,  there  will  be  a 
discussion  of  topological  considerations  in  relation  to 
the  classical  integral  theorems  of  vector  analysis*   These 
results  will  next  be  used  to  formulate  compatibility 
conditions  in  the  statement  of  existence  theorems  for 
harmonic  vectors  and  related  fields.   The  theory  of  these 
fields  will  be  pursued  further  to  reach  a  variational 
formulation  of  the  problems,  an  especially  valuable  approach 
for  fluid  magnetics.   We  present  also  a  general  existence 
theorem  which  goes   somewhat  beyond  our  immediate  needso 

V,      The  corresponding  Nevimann  problem  admits  only  constants 
as  solutions,  hence  B  is  the  gradient  of  a  constant. 
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We  present    also,   for   Its   own  intrinsic   interest,    a   treat- 
ment of   the   two-dimensional  case  which  is  not  merely  a 
specialization  of  the   three  dimensional  case,   but  owing 
to  the   difference   in  topology  displays   a  remarkable 
simplicity  and  elegance  of  its  own.     Employing  the  results 
for  haarmonic  fields  we   shall  exhibit   a  number  of  relevant 
well-posed  problems  for  the  pre-Maxwell   system  and  make 
some   application  of   the  methods    to  Faraday's  Law,   homo- 
polar  generators,    and   the  pinch   effect.      In  conclusion, :, 
there  is  a  brief  expository  comparison  between  the   Maxwell 
theory  and   the   theory   of  the  pre-Maxwell   equations. 
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(2.1) 


2.  Topological  Considerations 
The  classical  Integral  theorems  of  Gauss  and  Stokes 
form  the  most  natural  avenue  of  approach  to  the  formula- 
tion of  well-posed  boundary  value  problems  for  the  systems 
(l»6)-( 1.7)0   In  each  case  the  Integral  theorem  expresses 

the  value  of  an  Integral  on  a  manifold  In  terras  of  an 

7 
Integral  on  the  boundary,  viz,' 

a,  /  grad  0  *  dx  =  0{x^)   -  0ix.  ) 
C  "^      -^ 

b,  J     curl  Z  •  dS  =  /  Z  •  dx 

z:  r 

c ,  /  dlv  Z  dV  =  /  Z  •  dS 

D  S 

The  compatibility  conditions  mentioned  In  the  introduction 
are  restrictions  upon  the  botondary  values  consistent  with 
the  given  differential  equations  \inder  (2cl).   For  example, 
if  dlv  Y  =  0   is  one  of  the  differential  equations,  then 
Gauss'  formula  (2.1c)  implies  that  the  normal  component  on 
the  boxindary,   Y   ,   must  satisfy  the  condition   jp  Y  dS  =  0. 
In  the  same  way,  Stokes'  theorem  places  a  restriction  on 
the  tangential  boxindary  component  X.   of  a  vector   X 
which  satisfies   curl  X  =  0.   The  compatibility  require- 
ments can  be  formulated  precisely  only  in  terms  of 

TT   Here,  and  in  the  sequal,  we  systematically  use  the 

rbols,   G   for  an  open  arc,   P  for  a  closed  curve, 
for  an  open  surface,   S  for  a  closed  surface,  and 
D  for  a  three-dimensional  domain. 
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topological  properties  of  the  given  dDmainJ  specifically, 

in  terms  of  homology. 

Homology  of  curves  and  surfaces.   Relation  to  vacu-gm  fields. 

Let  us  take  a  three-dimensional  domain  D  and 
consider,  for  example,  the  application  of  Stokes'  Theorem 
(2.1b)  to  an  arbitrary  Irrotatlonal  vector  Z  in  Do 
If  the  surface  ^Z  which  is  bounded  by  the  closed  curve 
P  lies  entirely  within  the  domain  D,   then  /  Z  <>  dx  =  0. 
However,  a  closed  curve   ["*  may  exist  which  does  not 
bound  any  surface  JZ  Iji^S  entirely  in  D.   In  that 
event  we  may  not  conclude  that   ^  Z  «  dx  =0  for  every 
closed  curve   P  merely  because  Z  is  irrotationalo 
Nonetheless,  to  each   P  which  does  not  bound  in  the 
domain  D  we  may  associate  the  class  of  all  oriented 
closed  cxorves   P   such  that   P  -  P   does  bound  a 
surface  ^  in  D.   Consider,  for  example,   P  and   P 


Figure  1.   Closed  curve  homology  in  the  exterior  of  a 
torus. 
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in  Figure  1  where  D  is  the  exterior  of  a  torus.   It 

is  clear  for  all  curves  of  the  class  that 

J^ .    Z»dx=^  Z»dx  „  The  value  of  the  line  integral 

r  r 

therefore  does  not  depend  upon  the  curve  \    ^     but  only 

upon  this  class,  the  homology  class  of   f^^  Such  a 

constant  of  a  homology  class  is  called  a  period  of  the 

vector  field.  The  physical  meaning  of  these  results  is 

clear.  If  the  domain  D  is  interpreted  as  a  vacuum 

field  external  to  a  current  carrying  loop  then  the  m«m»f« 

(magnetomotive  force)  J^  B  -dx  around  a  closed  curve 

n  in  D  represents  the  current  in  the  loop  liriked 

by  the  given  cvirve.   This  current  is  a  period,  clearly 

independent  of  the  precise  choice  of  the  curve  and 

dependent  upon  homology  class  only;  in  other  words, 

the  ra«m»f .  is  the  same  for  all  closed  paths  which  link 

the  loop  once  in  the  same  direction. 

In  an  analogous  fashion  the  divergence  theorem  will 

yield  a  concept  of  homology  for  surfaces  when  applied  to 

solenoidal  vectors.   Consider,  for  example,  a  vacuum 

domain  D  between  two  concentric  spheres   (Figure  2), 

Two  surfaces   S   and  S   belong  to  the  same  homology 

class  if   S  -  S   boxinds  a  subdomain  of  D,   In  this 

case,  all  spherical  surfaces  concentric  with  the  inner 

sphere  belong  to  the  same  class.   Clearly  the  divergence 

theorem  implies  that   /^         /» 

/,  Z  *dS  =/  Z  -dS   , 
3         S 
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Fig  lire  2.   Closed  surface  homology  between  two  spheres. 

t 

where   S   and  S  lie  In  the  same  homology  class.   In 

this  example,  the  electric  flvLX   through  any  sphere  In  D 
containing  the  excluded  inner  sphere  is  a  period  and  is 
proportional  to  the  total  charge  exterior  to  D  within 
the  inner  sphere. 

Incompressible  fluid  flow  yields  another  example 
of  a  solenoidal  vector  field.   The  period  ^u  •  dS  ,   where 
u  is  the  velocity  vector,  is  assigned  to  the  homology 
class  of  a  closed  siirface  S   and  is  a  measure  of  the 
soiH'ces  contained  in  S. 

We  now  proceed  to  expand  somewhat  on  these  concepts. 
However,  so  brief  a  description  of  the  concept  of  homology 
as  the  present  one  must  necessarily  be  quite  incomplete. 
The  reader  is  referred  to  standard  works  on  topology  for 
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a  thorough  treatment  • 
The  homology  groups. 

The  three-dimensional  domain  D  is  assumed  to  be 
orlentable  and  to  possess  a  boiindary  consisting  of 
finitely  many  regular  bounded  closed  surfaces  •   In 
addition,  we  shall  assvime  regularity  on  the  part  of  all 
curves  mentioned.  For  example,  D  may  be  the  exterior 
or  the  interior  of  a  torus,  the  region  between  two 


Figure  3«   The  Borromean  rings. 


Wl      Alexandroff,  P.  S.,  and  Hopf,  H. ,  Topologie,  Berlin, 
Springer,  1935* 

Lefschetz,  S.,  Introduction  to  Topology,  Princeton 
University,  19i|.9. 

Seifert,  H.,  and  Threlfall,  W.,  Lehrbuch  der  Topologie, 
B.  G.  Teubner,  193il-» 

9,   A  domain  is  understood  to  be  an  open  connected  set.   For 
the  purpose  of  the  integral  theorems,  plecewise 
smoothness   is  sufficient  for  regularity.  For  topologi- 
cal purposes  it  is  required  that  the  surface  neighborhood 
of  a  point  on  the  boundary  is  topologie ally  equivalent  to 
a  disc.   If  the  domain  is  infinite  the  sphere  at  infinity 
is  to  be  thought  of  as  part  of  the  boundary,  and  in  that 
case  boundary  conditions  are  replaced  by  regularity 
conditions.   We  shall  not  generally  consider  other 
infinite  bounding  surfaces. 
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concentric   spheres,    the   exterior  of  the  Borromean  rings 
(Pigiire   3),   etc.      In  the    following,   we   assviine   a  finite 
topology,    explicitly  that  there   are   only  finitely  many 
independent  homology  classes  of   each  dimension.      To 
explain  this  concept,   let  us  recall   the    definition  of 
addition  of  curves   in  integration   theory  iiAiich  permits 
us  to   write 

/     -  I  -J 

^l^Z  ^1    °2 

(and  similarly  for  surfaces).   Each  curve  C  has  an 
additive  Inverse   -C   obtained  by  reversing  its  orienta- 
tion.  Clearly,  by  the  Introduction  of  a  zero  clement, 
the  set  of  curves  may  be  seen  to  form  an  additive  group 
which  contains  the  smaller  additive  group  of  closed  curves. 

A  closed  curve  is  said  to  be  homologous  to  zero  if  it 
bounds  a  surface  ^Z  ^^  ^*  '^^^   concept  of  bounding  is 


Figure  Ij..  Boundary  of  an  orientable  surface. 

10.  Device  on  the  coftt-of-arms  of  the  Milanese  family 
of  Borromeo  (114.06-),  thought  to  be  ancestors  of 
the  Ballentines, 
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roughly  defined  in  the  following  way;  A  closed  curve 
P  is  said  to  bound  a  surface  ^  if  it  constitutes 
the  entire  boundary  of  ^I  ^^'^  ^^   ^^   i^  possible  to 
partition  ^I  Into  non-overlapping  simple  surface  elements 
and  orient  their  boundaries  so  that  the  s\ira  of  the  bound- 
ing curves  of  the  surface  elements  la  the  curve   p. 
Thus,  we  see  in  Figure  ij.  that   P  -  P   bo\inds  T" 
while   P  +  P   does  not.  In  the  domain  D  of  Figure  1, 
the  curve   P  -  P   is  homologous  to  zero  and  we  write 
P  -  P  --^  0»  Two  cxirves   P  end  P   are  sal'  to  be 
homologous   ( P  '^  P  )  if  their  difference  bounds.  We 
may  then  speak  of  the  addition  of  homology  classes  rather 
than  curves,  where  the  additive  zero  is  simply  the  class 
of  all  curves  which  are  homologous  to  zero.   Integration 
over  a  linear  combination  of  cxirves  is  defined  by 

J  Z  •  <3lx  =  %  J     Z   o'dx  +  li.  /  Z  •  dxo 

%Q^-*V.C^  ^1  ^2 

If  curl  2  =  0,   then  J^    Z  '  dx     is  seen  by  this  formula 
to  be  a  linear  function  on  the  linear  vector  space  of 
homology  classes  of  closed  curves.   This  linear  vector 
space  is  called  the  homology  group  of  closed  curves. 
The  basis  of  the  space  is  restricted  to  be  finite.   In 

the  exterior  of  a  torus  (Figure  1),  the  curve   P 

11 
comprises  the  entire  basis r.  • 

11.  We  shall  use  any  of  the  members  of  a  homology  class 
to  designate  that  class. 

-  27  - 


The  homology  group  oj  slosed  siirfaces  is  defined 
in  an  entirely  analogous  ;  y»  A  surface  '   is  said  to 
be  homologous  to  zero  if  ,  t  bounds  a  subdc  lin  of  D» 
In  the  case  of  the  Interlc^  of  a  torus,  e\  y  closed 
surface  bo\inds  and  the  hoi.^logy  group  is  t  ;  null  space. 

Consider  the  example  of  the  Interior  *  .'  a  sphere 
with  the  interiors  of  two  coruses  deleted  (Figure  5)» 


Figure  5»   Homology  in  the  interior  of  a  sphere  with 
two  toruses  deleted* 


As  a  basis  for  the  closed  curve  group,  we  may  talce  the 

indicated  cijrves   1"^,   and   f^p.   For  the  closed  surface 

^    "it- 
group,  we  may  take  the  surfaces  S,  ,  Sp  of  the  toruses. 

One  further  point  of  which  we  shall  m&ice  signifi- 
cant use  for  the  purpose  of  assigning  compatible  boixndary 
values  in  the  differential  problems  (1.6),  (1.7  )J   it 
is  possible  to  take  curves  lying  on  the  boundary  as 
representatives  of  the  homology  classes.   In  Figure  5» 
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for  example,   we  may  replace      P^      and      P^     by      Pf 
and      Pg     on  the  boundary     S  ♦     We   remark  also   that   It   is 

possible    to   consider  the  homology  of  closed  ciirves  with 
respect   to   the    two-dimensional   domain  S    ;      i.e.,    the 
criterion  for  homology  is  whether   a  cxirve      P     bounds 
a  subdomain  of     S   •     We   observe  that   closed  curves  which 
do   not  bound  in     D     cannot  boxmd  on   the    two-dimensional 
boundary     S   .      However,    it   is   also   possible    to    take 
representatives   on     S       of   the  homology  group  of  closed 
curves   in  the    three-dimensional  exterior     D     of     D» 

In  fact,    it  can  be   proved   that   the  homology  group  of 

It 

closed  curves  on  S   is  the  direct  sum  of  the  homology 

groups  in  D  and  D. 

The  homology  group  of  closed  surfaces  is  more 
easily  described  than  the  group  of  closed  curves.   Let 
the  boundary  S   of  D  consist  of  q+1   closed  s\irfaces 
Sq,   S-,,..,,  S   ,   oriented  so  that   S-.  +  S-,  +  •••  +  S  --^0, 
the  sum  bounding  D»   Since  D  may  be  considered  as  the 
domain  common  to  the  interior  of  one  boundary  component 
and  the  exterior  of  the  others,  clearly  a  basis  for  the 
homology  group  can  be  foi^med  from  any  q  of  the 
components  of  the  boundary* 

It  Is  convenient  also  to  define  the  zero -dimensional 

12 •   Starred  manifolds  will,  without  exception,  refer 
to  manifolds  on  the  boxindary  of  the  domain  D. 
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homology.   Two  points  are  said  to  be  homologous  If  they 
bound  an  arc  in  D.  Homology  of  points  in  the  domain 
D  is  of  no  importance  since  D  is  connected  and  every 
pair  of  points  may  be  joined  by  an  arc.  The  zero- 
dimensional  homology  will  be  significant  here  only  for 
the  exterior  D  to  D,  which  has  as  many  separated 
components  as  the  common  boundary  S  . 
The  analytical  significance  of  homology* 

Let   Z  be  a  regular  vector  function  defined  in  D. 
To  each  closed  curve  and  each  closed  surface   3,   we 
assign  a  real  number  by 

[P]  =  [PjZ]  =  /z-  dx 

(2.2)  ,   ^ 

[S]  =  [S;Z]  =  j^  Z  «  dS   . 

If  Y  is  solenoidal,   div  Y  =  0,   then,   [S;Y] 
vanishes  for  any   S  which  is  homologous  to  zero. 
Consequently,  the  function   [SjY],   which  is  defined 
on  the  class  of  closed  surfaces,  has  the  same  value 
for  each  surface   S   in  a  given  homology  class.   The 
constant  associated  in  this  way  with  a  homology  class 
(rather  than  with  a  particular  surface)  will  be  called 
a  period.   Similarly,  if  X  is  irrotational,  curl  X  =  0 
in  D,   then  the  fiinction   [CjX]   is  defined  on  homology 
classes  independently  of  the  choice  of  individual  ctu*v#s. 
Conversely,   div  Y  =  0   (or  curl  X  =  0)   only  if 
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[SjY]   (or   [CjX] )   is  constant  on  each  homology  class. 
It  is  true  also  that  S,   and  Sp   (C,   and  C^)   are 
homologous  if  and  only  if   [S-|_}Y]  =  [S2>Y] 
(CG,jX]  =  [CpjX] )  for  every  solenoidal  Y  ( irrotational 
X).   It  would  have  been  possible  to  define  homology  of 
two  closed  curves  (surfaces)  by  requiring  that  the  two 
integrals   [  f^jX]   ([S}Y])   are  the  same  for  each 
irrotational  X   (solenoidal  Y). 

The  condition  curl  X  =  0  guarantees  the  existence 
of  a  single-valued  scalar  potential  locally.   In  order 
that  it  be  possible  to  define  a  potential  i)      over  all 
D  with  X  =  grad  j2$   it  is  both  necessary  and  sufficient 
that  the  closed  curve  periods   [ P;X]   all  vanish. 
Similarly,  for  any  solenoidal  vector  Y  a  vector 
potential  A   exists  locally  with  Y  =  curl  A.   The 
necessary  and  sufficient  condition  that   A  be  single 
valued  in  the  whole  domain  D  is  that  the  closed  surface 
periods   [SjY]   all  vanish. 

The  topology  has  been  assumed  finite}  that  is,  each 
closed  curve  or  surface  is  homologous  to  a  linear  combi- 
nation of  a  given  finite  set  of  basis  elements.   Let 
r^   ...  P'   denote  a  basis  of  the  closed  curve  group, 
S,  ...  S    a  basis  of  the  surface  group.   For  every 
closed  ciirve   f"'  and  closed  surface   S  we  have  the 
homology  representations 
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(2.3) 

11    c.  2.  q  q 

When  curl  X  =  0,   div  Y  =  0,   It  follows  that  all  periods 
are  defined  in  terms  of  the  finite  set  of  periods  of  the 
basis  elements  by 

[r;x3  =  a^^cr^jx]  +  ...  +  apCTpjx] 

(2.U) 

[s;y]   =bT[S.  JY]   +  ...    +b^CS^jY] 

XI  q     q 

It   is   Intuitively  clear  that   an     X     or     Y     can  always 
be  found  which  will   assign  arbitrary  given  values   to    the 
set  of  periods   of  a  basis.      Furthermore,    all   that  need 
be  known  about    any  solenoidal  vector     Y      (or  irrotational 
X)      in  order   to    determine      [S]      (or      [C])      for   any  closed 
surface    (or  curve)   are   the  values  of   its   periods   on  a 
basis. 

Consider   the   example    of    two   linked   toruses  bounded 
by   surfaces      S!'      and     Sp      as   indicated   in  Figure   6.      The 
domain     D     is   the   exterior   to   both.      The  oriented  closed 
curves      P,  ,    p,  ,    V^^*    l~2      ^^^  chosen   in  the    surfaces 
S-      and      Sp      and      P-      links  neither      S,      nor      Sg. 
Clearly      P_   /— '  0     while   no   one   of      P^,    P^,    p2»    '    £» 
separately,    bounds   in     D.      However      P^    — ■  \^ ^      and 
Pp   <^     P-      since   each   pair  boxinds   an  annular  disc. 
Clearly      P,      and      \~^^      could  be   taken  as   a  basis   of  the 
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Figure  6 «   Homology  In  the  exterior  of  two  linked  toruses, 

closed  curve  group  in  D«   There  exists  a  two -parameter 
infinity  of  homology  classes  constructed  so  that  every 
closed  curve  in  D  is  homologous  to  rar,  +  npV   for 
some  m  and  n.   In  this  example,  the  closed  siafface 
group  has  a  basis  consisting  of  two  elements  which  may 
be  taken  as   S,   and   Sp.   (The  third  component  of  the 
bo-undary  is  the  sphere  at  infinity).   Note  that  homology 
alone  does  not  distinguish  between  this  case  and  that 
of  Figure  5  although  the  two  are  topologically  distinct. 
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3»   Homology  Modulo  the  Boundary* 

The  homology  of  closed  manifolds  has  just  been  seen 

to  be  intimately  connected  with  certain  differential 

operators <  the  closed  curve  homology  with  the  curl,  the 

closed  surface  homology  with  the  divergence.   For  an 

irrotational  vector  X  we  have  observed  that  the  integral 

about  a  closed  curve   P,   [  PjX]  =  ^  X  •  dx,   depends 

only  on  the  homology  class  of   P;   for  a  solenoldal  Y, 

the  integral  over  a  closed  surface  S,   [SjY]  -  v   Y  t  dS, 

S 
again  depends  only  on  homology  Class,  We  shall  find  that 

when  the  homogeneous  differential  relation  is  supplemented 
with  a  suitable  homogeneous  boundary  condition  the 
relevant  periods  all  vanish  but  other  integrals  arise 
which  again  depend  only  on  the  topology  and  not  on  the 
particular  choice  of  the  manifold  of  integration. 

Consider,  for  example,  a  domain  D  and  an  irrota- 
tional vector  X  in  D  with   X.  =  0  prescribed  on 
the  boundary   S  .   Since  a  basis  for  the  homology  group 
of  closed  curves  in  D  can  be  chosen  on  S  ,   it  is 
clear  at  once  that  [ P}X]  =  0  for  all  closed  curves 
in  D.  Consider,  now,  two  arcs  C   and  C   which  begin 

and  terminate  on  S  •   If  it  is  possible  to  find  arcs  on 

t 
the  boundary  which  taken  toget^ier  with  C   and  C   form 

a  closed  curve   P,   then  clearly,  under  the  assumption 

X.  =  0,   the  only  contributions  made  to   [  P;X]   come 
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from  C  and  C  ,   It  follows  that  the  integrals  on  C 
and  C   are  equal j  the  open  arcs  C   which  with  C 
can  be  completed  to  a  closed  curve  by  additions  on  the 
boxindary  constitute  a  class  upon  which  the  line  integral 
J  X  •  dx  is  constant* 

c 

This  suggests  a  second  kind  of  homology  which  we 
shall  call  homology  modulo  the  boundary.  This  new  concept 
of  homology  is  most  clearly  understood  by  relating  the 
homology  of  the  domain  D  to  that  of  the  exterior  domain 
D  complementary  to   D. 
The  Alexander  duality  theorem. 

There  is  a  duality  relation  between  homology  in  D 
and  homology  in  the  exterior  domain  D  which  is  given 
by  a  remarkable  general  theorem  of  Alexander.   For  a 
three  dimensional  space,  in  particular,  the  Alexander 
duality  theorem  states  that  the  number  of  independeait 
elements  of  the  k-dimensional  homology  group  of  D  is 
equal  to  that  of  the   (3  -  k  -  1 )-dlmensional  group  of  D. 
The  two  groups  are  called  dual  to  each  other. 

Prom  the  duality  theorem,  it  follows  that  the  closed 
curve  group  in  D  is  dual  to  the  closed  curve  group  in 
D.   As  we  have  remarked,  the  homology  group  of  closed 
CTxrves  on  the  common  boundary  S   is  the  direct  sum  of 
the  groups  in  D   and  the  exterior  D.   We  conclude  that 
the  number  of  linearly  independent  closed  curves  on   S 
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.Po 


Flgore  7«   Dual  bases  in  a  sphere  with  two  toruses  deleted. 


is  double  the  number  of  those  in  D.  For  example,  if  D 

is  the  interior  of  a  sphere  with  two  toruses  deleted 

depicted  in  Figure  7,  V^-^      and   P^  constitute  a  basis 

for  D,   and   P^   and   P^  ^°^  ^«  "^^^  basis  of  the 

■if- 

closed  curve  group  on  S   consists  of  representatives 
from  each  of  the  four  homology  classes. 
Homology  of  open  curves  modulo  the  boundary. 

The  dual  to  the  homology  of  closed  svirfaces  is 
easily  visualized.  Let  the  bovmdary  S  of  D  be 
written  as  the  sxim  of  its  components. 


a:    &       + 


S^  -*-•••  -t-S,   To  each  component  S,   we  may 
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assign  a  point     P^     in  the  component  of  the  exterior 
D     bounded  by     Sj   (Figure    7).     The   zero -dimensional 
homologn-   group  of     D     is  generated  by   the  basis      P^P,,* 
(k  =  l,2,...,q). 

If     P.^     and     P.      are    taken  on   the   boundary  they 
it-aj  be   connected  by   an   arc      C,       in     D,      (Figure   8), 


Fig'-ire   S«      Representation  of   the   open  arc  homology. 


The  arcs  C,  may  be  tal'ien  as  a  basis  for  the  homology 
in  D  of  open  curves  r.odulo  the  bo-'-Lndary.  Two  curves 
C      and      C        in     D     which  initiate   and  terminate   on     S 


are    said  to   be  hornoloscus  nodulo   the   bo'jLndary 
(C    ■ 


1  jft. 

C         [mod  S    J )      if  their  zero-dimensional  botmdaries 


are  homologous   in     D.      It   is    completely  equivalent   to 
set      G    ~    C        [mod  S   ]      when   it   is  possible   to    complete 
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C    -  C        to   a  closed   curve  which  bounds   in     D     by  adding 

arcs   on  the   boundary  only.      The  homology  group  of  open 

arcs  modulo   the  boimdary  in     D     is   directly  isomorphic 

to   the  zero -dimensional  homology  group   in  the   exterior 

D.     Alexander's  duality  theorem  therefore   enables  us  to 

establish   a  duality  between  two   homology  groups   in     D, 

the   group  of  closed  surfaces     fSj,      and   the  group  of 

open  ciorves      |CJ» 

The  bases      3,       and     C,  ,      (k  =  1,    ...,    q),      of  the 

dual   groups   chosen  above    are  reciprocal ;    i.eo,    the  number 

of   intersections    (an  intuitive   concept  whi'jh  can  be 

defined  precisely)   of     C.      and     S,       is   equal   ■'"o     1     if 

J        '^ 

j  =  k   and  equal  to   0   otherwise o 

Homology  of  open  surfaces  modulo  the  boundary* 

The  duality  between  the  bases  of  the  closed  curve 
groups  in  D  and  D  is  somewhat  more  complex.   It  can 
be  shown  that  there  are  dual  bases,   r,   in  the  domain 
D,  and   n*.  in  the  exterior  D,   ( j,k  =  1,  2,  .» • ,  p) 
with  the  property  that   P,   links   ["*,   exactly  once 
if   J  =  k  and  not  at  all  if   j  /^  k   (cf.  Figure  7). 
Moreover,  by  taking  the  representative  curves   ("I 
and  \.      on  the  boundary,   S  ,   we  can  obtain  in  D 
a  duality  between  the  closed  curve  homology  and  the 

homology  of  open  surfaces  modulo  the  boundaJ*yo   The  curves 

A  -  4 

\,       do  not  bound  in  D  or  on  S  ,  but  can  be  chosen  to 
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bound  surfaces  y,       in  Do  Each  surface  "T"^  will 
have  the  property  that  the  intersection  number  of  the 
closed  curve   P.  with  the  open  surface  7~t.  is  1 
if   j  =  k  and  0  otherwise*   The  bases   p.,  T", 

J 

chosen  in  this  way  will  then  be  reciprocal*  The 
homology  of  the  open  surfaces  y,  is  defined  above  by 
the  homology  of  the  boundary  closed  curves  P,  but  it 
is  intuitively  clear  that  this  is  equivalent  to  setting 
y~  "  y  '^  0  [mod  S  ]  if  it  is  possible  to  complete 
YZir  "  IZir  ^°  ^  bounding  closed  surface  in  D  by  adding 
a  part  of  S  « 

An  instructive  example  is  that  of  the  exterior 
domain  to  two  linked  toruses  (Figure  6)o   The  curves 
p.  ,  Pp   serve  as  a  basis  of  the  closed  curve  homology 
in  D  while   P,   and  \~^^      serve  as  a  dual  basis  in  D« 

A.  A 

The  closed  curves   P,   and  [''^      do  not  bound  in  D, 

y\  A 

but  the  curves  ^±   ~    ^2    »       ^Z   ~    ^1     ^"  their  respec- 
tive homology  classes  boimd  annular  domains j,   >  ..  p  ^p* 
The  open  surfaces  ">"■.  ,  y^   constitute  a  basis  for 
the  homology   [mod  S  J   which  is  reciprocal  to  the  basis 

'1  »  '2° 

The  transformation  rule  for  reciprocal  bases  is 

easily  obtainedo  If  P^  and  >^^,  (j»k  =  l,2,ooo,p), 
are  reciprocal  bases  and  \.  is  a  second  basis  for  the 
closed  curve  group  in  D,   then  a  reciprocal  basis 
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«  t 

/y.   9      to   I  <  >   Is  obtained  from  the  J~,   by  the  inverse 

transpose  of  the  linear  transformation  which  takes  the 

n*.  into  the   p..   We  shall  restrict  o\arselves  to 

reciprocal  bases  henceforth. 

The  analytical  significance  of  homology  modulo  the  boundary^ 

In  the  case  of  a  solenoidal  vector  Y,   ws  know 

that  the  surface  integral  ^  Y     dS  is  defined  as  a 

S  ^ 
function  of  homology  class  alone  and  does  not  depend 

upon  the  special  choice  of  the  surface  S.   If  the  differ- 
ential relation   div  Y  =  0   is  supplemented  by  the 
condition  that  the  normal  component  of  Y  vanish  on  th© 
boundary,   ^n  ~  ^   °^  ^  »   then  the  periods   [SjY] 
vanish  on  the  basis  S^,   (k  =  1,2, ...,q),   chosen  on 
S   and  consequently  all  closed  surface  periods  vanish* 
In  that  event  it  becomes  possible  to  define  open  surface 
periods  in  terms  of  the  homology   [mod  S  ]•   If  ^  is 
a  surface  in  D  with  boundary  on  S   we  set 

(3.1)  CII;Y3  =  f  Y.  dS  . 


The   fact    that   the   period      [^IjY]      is    a  function  of 
homology   class    alone   is    clear    if  we  observe   for   any 
T~     --^  ^^Z.     [i^od  S   ]      that     ^  -  ^        can  be   completed 
to    a  closed  boxonding   siorface  by  the    addition  of  part 
of     S   .      The   converse   of  the   stated  result   is  also    true. 
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If  the  Integrals   /  Y  •  dS  are  periods,  I.e.,  If  they 
depend  upon  homology  class  alone,  then  the  closed  surface 
periods   [S;Y]   vanish  and  Y  =  0  on  S' . 

These  properties  can  also  be  stated  in  terms  of  the 
vector  potential.   Since  div  Y  =  0   there  exists  a 
vector  potential  A,   defined  locally,  for  which 
Y  =  curl  A.   If  the  open  surface  integrals  are  periods, 
there  exists  a  single -valued  vector  potential  A,   where 
A  is  a  local  surface  gradient  on  S    .  This  result 
is  plain  from 

A  A  •  dx  =  S^   curl  A  •  dS  =  A.  Y„  dS  =  0   , 

r*         r:^  ^  n 

where      >  is    any  surface  element   on     S'       and      \^     is 

its  boTxndary.      If  it   should  happen  further  that   all 
periods    V^ZX^^      vanish  then     A      Is   the    surface   gradient 
of    a  single-valued  potential   on     S   . 

Similar  results    to    the   above   can  be   obtained  for 
an  Irrotatlonal    vector     X.      Prom     curl  X  =  0      it   follows 
that   the   closed  curve    Integral      ^^  X  •  dx      is   a  period. 

r 

If  the  boundary  condition  is  that  the  tangential  component 

vanish,   X.  =  0   on  S  ,   then  the  periods   [  njXJ 

13»   For  the  purposes  of  this  paper,  a  vector  A  is 

defined  as  a  local  siirface  gradient  on  the  sizrface 
S   if  J^A»dx=0  for  any  closed  curve   f^ 

r 

bounding  on     S.      We  write      A  =    '^«0     where     ^(x)    = 

X  ^ 

r  A*  dx,   the  path  of  integration  being  taken  in 

the  surface, 
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vanish  on  the  basis   f"^  »   (k  =  l,2,...,p),   chosen  on 
S   and  consequently  all  closed  cxirve  periods  vanish* 

The  open  curve  integrals  for  curves  C  with  boundary 

It 

©n  S   now  become  periods, 

(3.2)  [CjX]  =  J  X  •dx     . 

C 

Conversely,  if  the  integrals  /  X  •  dx  are  periods  then 

C 
the  closed  curve  periods   [  f^jX]   vanish  and  X.  =  0  on 

Since  curl  X  =  0  there  exists  a  local  scalar 
potential  0     with  X  =  grad  0*      If  the  open  curve 
integrals  are  periods  then  0     is  single-value i  in  D 
and  constant  on  each  coinponent  of  S  •   If,  In  addition, 
all  values   [C|X]   vanish,  then  0      takes  the  same  constant 
value  on  every  component  of  S  • 

In  each  case  there  actually  exists  a  vector  which 
realizes  arbitrarily  specified  period  conditions  on  a 
basis.   For  example,  if  D  is  a  vacuum  domain,   div  E  =  0. 
We  may  freely  specify  the  charge  Q^^  =  lc[S,  ,B3, 
(k  =  l,2,..«,q),   in  the  exterior  domain  to  each  of  the 
elements  of  the  basis  of  closed  surfaces  taken  as 
components  of  the  boundary  with  no  restriction,  but  not* 
that  ^  Q,  =  0,   Since   div  B  =  0,   the  non-existence 
of  magnetic  charge  implies  that  the  periods   [SjB] 
always  vanish.   If  we  impose  the  condition  B  =  0  on 
S   appropriate  to  a  perfectly  conducting  boundary  we 
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observe  that  the  open~surface  periods  are  Interpretable 
as  magnetic  fliuces.   These  fluxes  X  =  E^,  SB] , 
(k  =  l,2,o»»,p)»  may  be  given  arbitrarily  by  specify- 
ing  ctir rents  in  external  circuits   p,   lying  in  the 
homology  classes  of  the  boundaries  of  the  respective 

For  the  purpose  of  convenient  reference,  the 
results  concerning  the  relations  between  analytical 
and  topological  concepts  are  summarized  in  the  follow- 
ing table »   Periods  are  said  to  be  veil-defined  if  the 
integrals  are  determined  by  homology  class  alone© 
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Table  I 
Summary  of  the  Relations  of 
Analytical  Concepts  to  Homology 


la 
i 

ii 
ill 

iv 


curl  X  =  0 

[  P;X]   well-defined 

[  TiX]  =  0 

X  =  grad  f) 

curl  X  =  0    , 
^t   =° 

[CJX]   well-defined 

X  =  grad  ^   , 
)Zf  =  )Zi^  on  Sj; 

[CJX]    =  0 

X  =  grad  ^   , 
JD  =  const,   on  S 

lb 
i 

ii 
iii 


iv 


div  Y  =  0 

[SJY]   well-defined 

• 

[SJY]    =  0 

Y  =  curl  A 

div  Y  =  0    , 

[2I»Y]   well-defined 

Y  =  curl  A   J   A  is 

Y^=0 

local   surface 
gradient   on  S 

[JIjy]  =  0 

Y  =  curl   A   J 

^   single -valued 
on  S 

The  table  is  to  be  read  as  follows.   In  each  row  of  the 
table,  the  separate  entries  are  equivalent.   Row  i  in  la, 
for  example,  should  be  read,  "Closed  ciorve  periods   [  PSX] 
are  defined  as  functions  of  homology  class  alone,  if  and 
only  if  curl  X  =  0  in  D."   As  one  reads  down  a  column  of 
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either  table  the  conditions  become  Increasingly  more 
stringent  so  that  each  of  the  entries  in  a  row  implies 
all  the  entries  on  the  rows  above?   so,  for  example, 
from  X  =  grad  f6   in  la,  ii  we  may  conclude  that 
curl  X  =  0  and  [  flX]     well-defined. 
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I|.»   Compatibility  Conditions* 
The  topological  concepts  developed  in  the  two 
sections  preceding  will  be  reflected  directly  in  the 
formulation  of  the  boundary  value  problems  of  direct 
Interest  to  us.  Prom  the  analytical  properties  of 
homology  we  shall  derive  certain  evident  compatibility 
conditions  which  are  to  be  employed  in  the  formulation 
of  well-posed  problems  for  the  pre-.Maxwell  equations* 
Compatibility  for  values  of  the  tangential  ooxindary 
component  of  an  irrotatlonal  vector a 

Let  us  first  consider  an  irrotational  vector  X, 
curl  X  =  0   in  the  domain  D   and  see  what  restrictions 
are  imposed  on  the 'assignment  of  tangential  boundary 
values   X.   on  s'"'*   For  any  closed  cxirve   P^  which 
bounds  a  surface  element  2_"   °^   ^  »   ^®  must  have 
<r_^   X.  •  dx  =  Jj.   curl  X  •  dS  J   hence, 

(U.l)        4^  ^t  •  ^  "^  °  (P^OonS*). 

If  the  line  integral  of  the  surface  vector  X^  vanishes 
for  every  curve   f""^  which  boxmds  on  S  ,   then  X  is  a 
local  surface  gradient  (cf»  Footnote,  p.  h,\)o      On  the 
boundary   s"*,   the  homology  group  of  closed  curves  has 
as  a  basis  the   2p  curves  X\    •*•  '  p*  '^1  *"  '  p  * 
where  V^    .,,    \^^     is  a  basis  for  the  homology  group 

1  P  A,  /\^ 

Of  closed  curves   in     D,      and      f^    -  •    F^     ^°^^   surfaces 


5~-i  o»>»  T"-  In  D,  which  constitute  a  reciprocal  basis 
for  the  homology  group  of  open  s\irfaces.   It  follows  that 

[  r^,X]  =^  X»dx  =  ^  X^«dx=y   curl  X  •  dS  =  0  • 


'  k 


r? 


We 


see,  then,  that  boundary  values  X.   cannot  be  assigned 


freely  to  an  irrotational  X  in  D,   but  that  half  the 
closed  curve  periods  of  X^^  on  the  boundary  S   must 
vanish? 


ik.Z)     [r^,X]  =  /  X^  -dx  =  0 

1   V 


(k  =  1, •••,?)  • 


For  an  irrotational  X  we  conclude  that  to  prescribe  the 

■it 

tangential  field     X^      compatibly  on     S   ,      the   tangential 

field  must  be   a  loftal    surface   gradient   and  half  its 
periods  must  vanish  while   the  remainder  are   completely 
determined  by   the   given  boundary  values. 

As   an  example   of  these   restrictions,    consider   a  helical 


Figxjre  9«     Helical  vector  field  on  the   surface   of  a  torus, 
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field  X^  on  the  surface  S   of  a  torus  (Figure  9)» 

Suppose  X.   to  be  a  local  surface  gradient.  There  are 

■Si- 
two  fundamental  closed  curves  on  S   indicated  in  the 

/,         /\ 
figure  by   P  and   P  where   P  boimds  in  the  interior 

domain  and   P  boimds  in  the  exterior.   To  be  compatible 

with  irrotational  X  inside  the  torus,  the  transverse  period 

[  P,X]  must  Vanish  (the  pitch  of  the  helix  Is  said  to 

be  infinite).  To  be  compatible  with  irrotational  X  • 

in  the  exterior  of  the  torus,  the  longitudinal  period 

[  P,X]   must  vanish  (the  pitch  of  the  helix  is  said  to 

be  zero). 

Next  consider  an  irrotational  vector  satisfying 

the  homogeneous  boxmdary  condition  X.  =  0.   All  closed 

curve  periods  vanish  and  the  compatibility  conditions 

are  trivially  satisfied.  This  time  it  is  the  open  curves 

which  become  of  interest.   Open  curve  periods  are  defined 

by  setting 

(i|.3)        [C;X]  =/  X^dx 

C 

where   C   is  an  arc  with  endpoints  on  the  boundary  S  • 
In  that  case,   X  =  grad  0,      where  ^   -  ^]^»      ^   constant, 
on  the  boiondary  component  S,  .   The  periods  on  a  basis 
are  the  potential  differences 

where   C,   is  an  arc  connecting  ST   to   Sq,   (k  =  l,2,.».,q). 
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The  conditions  cited  are  Indeed  compatible  In  the 
sense  that  there  always  exists  at  least  one  vector  X 
having  boundary  values  X^^  specified  arbitrarily 
subject  to  the  conditions  (i^.*!)  -  (ij..2).  Moreover,  In 
case  X.  =  0,  there  exist  vectors  X  for  each  possible 
specif ication  of  the  constant  values  of  the  potentials 
on  the  boundary  surfaces  S,  ,   (k  =  0,1,2,  ••.  ,q)» 
Coropt ability  for  values  of  the  normal  boimdary  compenent 

of  a  solenoidal  vector^ 

For  a  solenoidal  vector  Y,   div  Y  =  0  in  D,   the 

specification  of  compatible  normal  boundary  values  Y 

is  simple.   Let  the  boundary  S   be  written  again  in 

the  form  S  =  3-^  +  S,  +  •••  +  S   where  the  terms  are 

u    1         q 

the  separate  components  of  the  boundary.   Since  S  -"^  0 
in  D,   values  of  Y   will  be  assigned  compatibly  only  if 

(4.5)     ^  Y  dS  =^  [S*iY]  =  0   . 
S*  n      k=0   ^ 

Under  the  homogeneous  bo-undary  condition  Y  =0, 
all  closed  surface  periods  of  the  solenoidal  field  vanish, 
but  now  the  open  surface  periods  are  defined  in  terms  of 
homology  modulo  the  boundary.   If  Yl     ^^   ^   surface  in  D 
with  boxmdary  on  S  ,   we  set 

(U.6)  CIIJY]  =y  Y«  dS   . 

As  In  the  preceding  case,  the  conditions  are  compatible 
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in  the  sense  that  there  always  exists  a  aolenoldal  Y 
for  every  specification  of  Y   subject  to  Ik^S]*      In 
addition,  in  the  case  Y  =  0,   the  open  surface  periods 
(Ij.»6)  may  freely  be  assigned  on  a  basis  y"-  , «« .  -T"-* 
Compatibility  conditions  for  the  boundary  value  problems 
with  InhomoKeneous  differential  relations*    " 

If  the  differential  equation  is  Inhomogeneous, 
ciirl  X  =  ^     (with  dlv  ^  =0)  or  div  Y  =  (T,   then 
the  concept  of  period  as  defined  in  terms  of  homology 
classes  must  be  modified  and  the  associated  compatibility 
conditions  altered  accordingly* 

In  the  case  of  the  equation 

(i|.7)  div  Y  =  ^ 

where      C       is   a  prescribed  scalar  fxinction   in     D     we  have 
Instead  of    (1+.5)   the  compatibility  condition. 


(I|..8)  jZf    Y^  dS  =  /  0- 


dV 


S  D 


Furthermore,    the   Integral    J^  Y  •  dS     is  not   constant   on 

S 

the  homology  class  of   S,   Nonetheless,  if  S  ' — '  S   it 

is  possible  to  evaluate  the  si;irface  Integral  on  S   in 

terms  of  that  on  S  and  the  values  of  ^  ;   for  letting 

t  t 

D   denote  the  domain  bounded  by  S  -  S  we  have 


y  Y  •  dS  -  y  Y  •  dS  =  /  0'  dV  . 

s'       s      d' 
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Therefo2?e  we  set   [SjY]  =  _/ Y  *  dS  and  call  the  integral 

S 
a  period,  but  it  must  be  kept  in  mind  that  the  values  of 

the  integral  have  to  be  calculated  from  its  value  on  a 

specified  member  of  each  homology  class  and  from  the 

known  values  of  (T  •   In  prescribing  periods  for  an 

inhomogeneous  differential  system  we  must  exert  care  to 

take  fixed  representatives  in  each  homology  class. 

In  exactly  the  same  way  we  consider  the  equation 

(1|..9)  curl  X  =  ;^  ,  (div  ;=]  =  0)  , 

where  "^     is  a  prescribed  vector  in  D»   The  values  of 
jT  X  •  dx  do  not  depend  upon  homology  class  alone  but 
depend  also  upon  the  values  of  _^  •   Thus,  if   |^  ^^   |^ 
then   r  -  P  bounds  a  surface  ^  in  D  and 


^    X'dx-  /x»dx=/'  j::,«dS 


The  definition  of  closed  curve  period  must  therefore  be 
extended  in  terms  of  the  values  of  ^  •   Again,  specific 
values  of  the  periods   [ PjX]  =  $_X  •  dx     must  be  assigned 
to  fixed  representatives  of  the  homology  classes,  but 
when  one  such  value  is  assigned  all  others  in  the  same 
homology  class  are  also  determined. 

The  compatibility  conditions  on  the  boundary  values 
X.   are  also  altered.   No  longer  must   X.   be  a  local 
surface  gradient,  but  must  satisfy  instead  of  (ij-.l), 
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where  J~^'      is  an  arbitrary  surface  domain  on  S   with 
boundary  l' •      This  relation  amounts  roughly  to  the 
specif loatlon  of  the  surface  curl  of  X^  as  a  given 
function  on  the  surface j  It  determines  X^     uniquely 

to  within  an  arbitrary  additive  surface  gradient, 

-it 

Furthermore,  half  of  the  closed  curve  periods  on  S 

are  determined  by  the  values  of  _^  ;   thus,  Instead 
of  (I|.«2)  we  have 


•dS 


surface  y~,   bounded  by   PT.   since   dlv  _^  =  0.  )   These 


The  period   [  f~7  jX]   will  now  depend  upon  the  choice 
of  the  representative  of  the  homology  class.   (However, 
the  value  of  the  period  is  independent  of  the  particular 

-k  --— °-  -^   r^k 

rather  complicated  compatibility  conditions  may  be  described 
conveniently  in  the  following  way.   Let   X   be  any  vector 

which  satisfies   curl  X  =  ]^  •   The  vector  X^  has 

0       * 

certain  bovindary  values  X,   on  S   and  certain  periods 

(^^   X,  •  dx   on  the  specified  curves   pf  •   For  boundary 

■k  

values     X.       to   be   compatible  with  the   given  values     j^    , 

X.      must  differ  from     X,      by   a  local    surface   gradient 

and     X,      must   have   periods    equal   to    those  of     X.       on   the 
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The  concept  of  open  period  for  Inhomogeneous  boundary 
conditions. 

The  same  type  of  generalization  may  be  used  to  define 
periods  on  open  curves  and  surfaces  when  the  specified 
boundary  values  are  inhomogeneous.   For  example,  consider 
curl  X  =  0,  X.   given;  although  /  X^.  •  dx  is  not  a 
function  of  homology  class   [mod  S*]   alone,  the  values 
of  the  integral  for  two  curves  in  the  same  homology  class 
differ  by  the  value  of  a  curvilinear  integral  on  S 
which  is  known  from  the  values  of  X. .   Similarly,  if 
div  Y  =  0,  y   given,  the  surface  integrals  en  two  open 
surfaces  in  the  same  homology  class  differ  by  an  integral 
over  a  surface  domain  on  S   which  is  known  in  terms  of 
the  given  values  Y   • 

If  the  differential  relations  are  inhomogeneous 
also,  the  difference  between  the  open  periods  of  two 
elements  in  the  same  homology  class  will  have  an 
additional  term  defined  as  above  by  the  given  data  in  Do 
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5«  Existence  Theorems  for  Harmonic  Vectors 
and  Related  Fields. 

In  this  section  we  formulate  a  nxamber  of  well-posed 
problems  for  harmonic  vectors  and  the  related  vectors 
satisfying  inhoraogeneous  differential  equations.   The 
topology  of  the  domain  will  be  reflected  in  the  existence 
theorems  by  the  period  conditions  and  restrictions  on 
the  boundary  conditions* 
The  linear  decomposition  of  solutions. 

The  solution  of  the  pre-Maxwell  equations  can  be 
reduced  as  we  have  indicated  in  the  introduction  to  the 
successive  solution  of  systems  of  the  kind 

(5.1)         _  ^ 

{div  j^  =  0)    • 

It   is  convenient   to   consider   a   solution  of   system    (5*1 ) 
as   the   sum  of   solutions  of   the    systems 


(5.2) 


(5.3) 


(5.14.) 


(div    -    =  0) 


-  Bk  - 


The  system  (5«3)  has  the  particular  solution 
Z  =  grad  ^ 

D 
and  (5»ij-)  has  the  particular  solution 

Z  =  curl  A   , 

(5.6)         1  ri  — 

D 

If  H  represents  a  harmonic  vector,  that  is,  a 
solution  of  (5»2),"we  see  that  a  solution  of  (5.1)  may 
be  put  in  the  form 

(5.7)       Z  =  H  +  grad  fi   +  curl  A  • 

Since  0   and  A  are  explicit,  the  role  of  properly 
posed  boiindary  conditions  is  to  fix  H  uniquelyo 

On  the  other  hand,  any  given  dif ferentlable  vector 
Z  may  be  expressed  in  the  form  of  the  decomposition 
(5.7)  by  taking  for  cT   and  ^  in  {S^,^)   and  (5.6), 
div  Z   and  curl  Z,   respectively*   We  conclude  that 
any  vector  Z   can  be  written  as  the  sum  of  a  harmonic 
vector,  a  gradient,  and  a  curl,   A  slightly  different 
formulation  is  that  every  vector  can  be  written  as  the 
sum  of  an  Irrotational  and  a  solenoidal  vector   (H   is 
both).   In  general,  neither  the  irrotational  nor  the 
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solenoidal  vector  has  a  single-valued  potential;  the 
special  feature  of  (5.7)  is  that  the  multi-valuedness, 
if  any,  is  concentrated  in  H.   It  is  possible  to  refine 
such  decompositions  in  such  a  way  as  to  include  in  the 
statements  themselves  existence  and  uniqueness  theorems 
of  the  systems  (5.1)  -  {B»k)»      This  will  be  done  in  the 
next  section.   But  we  shall  first  develop  such  existence 
theorems  directly  in  the  language  of  partial  differential 
equations  and  in  the  next  section  restate  them  in  terms 
of  projections  in  a  Hilbert  space,  a  form  particularly 
appropriate  for  variational  analysis. 
Existence  theorems  for  harmonic  vectors. 

We  first  consider  harmonic  vectors,  i.e.j  solutions 
Z,   of  {B»2)    in  a  given  domain  D.   The  subtleties  are 
purely  topological.   In  a  simply  connected  domain,  we 
have   Z  =V^     where  A0  =  0  •   The  problem  for  harmonic 
Z   reduces  to  the  scalar  potential  equation  which  can 
be  solved,  for  example,  if  0     or  ^0/<^n      (subject  to 


conditions  on   Z  we  may  prescribe  either  the  tangential 


.  ,   dS  =  0)   is  given  on  the  boundary.   For  boundary 


component   Z.   as  an  arbitrary  surface  gradient  or  the 
normal  component   Z        subject  only  to  ^Z  dS  =  0  • 
In  a  more  general  domain  the  problem  becomes  one  of 
specifying  periods  as  well  as  boundary  data,  and  verify- 
ing their  compatibility.   Prom  the  results  of  the  last 
section,  a  harmonic   Z  always  has  well-defined  periods 
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^Z  •  dJC  =  [  p]   as  well  as   /z  •  dS  =  [S].  With  the 
boundary  condition  Z  =  0  ,   the  periods  l^]   =  Jz  •  dS 
become  meaningful  while  the  periods   [3]  vanish.   If 
Z^  =  0  ,  then  the  periods  [0]  =  J Z  •  dx     become  mean- 
ingful while  the  periods   [ P]  vanish.  Note  that  specific 

periods  [  p]   and  [T~]   (both  defined  when  Z„  =  0 )  are 

n 

not.    In  general,    equal   in  value   even   though  there  may  be 
a  correspondence  between   them  in   an  isomorphism  of   the 
groups     [Tl]     and     JP/. 

We  have  found  four  homology  groups    in     D,      closed 
surface   denoted  by     {St,      open  cxirve  by     |cf»      closed 
curve      l.n'j*      fiJ^d  open  surface     iZI|»     ^J  choosing 
reciprocal   bases  we  have  exhibited   a  duality  between 
the  closed   cxorve    group    j  Pj      and  the   open  s\irface   group 
j^3j»      ^^'^    another  between  the    closed   surface  group     ^SJ 
and   the   open  curve  group     {^j    •      I^^  terras   of   linear   spaces 
there    are    then  only  two    abstract   groups «      Of  the    two 
abstract  groups,    the      G,3     group   Is    the  more  often 
associated  with  the   electric   vector  and   will  be   denoted 
by      C-»      The   periods   of   the        L-    group   in  Its   two 
representations,    open   and   closed,   may  be   thought   of   as 
e.m.f.'s,      V  =   [CjE],      or  charges,      Q,  =  kCS;E]    o      The 
P,^     group   Is   usually  more    appropriate   to    the   magnetic 
vector   and    Is   denoted  by     ^^»      The   open    and   closed 
T^t-perlods    are    typically  defined  by   fluxes     J  =   [^TJ^^* 
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and  cur3?ents     jxl  =  [  n;B]      In  external  circuits.     We 
shall  on  occasion  refer  to   a   set  of  periods  of  these 
abstract  groups  as      Mj-periods  or      c -periods  without 
specifying  the  particular  representation  of   the    abstract 
group. 

The  general  well-posed  problem  will  give  boundary 
values   of   either     2        and  the   values   of   a  set   of     A>t -periods 
or     Z.       and  a   set  of       C-periods,     We    shall   see   that   the 
values   of   dual  periods      V^]      and      [  P*] ,      for  example, 
while   not   equal    are  related  in  the    sense   that   specification 
of  one,    together  with  boundary  conditions,    determines   the 
other   and   therefore   only  one    set  of  periods   i'rv   the      '^ 
group  may  be    fixed   independently.      A   similar   statement 
holds   for   the       C'     group. 

Without  proof,   we    state  ^J 
Theorem  1?      There   exists   a  unique  harmonic   vector     Z 
with  arbitrarily  assigned     /?t-periods    (in  either   the 
)y^ \     or      )\~^(      representation)    and  with    arbitrarily 
assigned  values     Z        on  the  boundary     S        subject   only 
to  the  compatibility  condition    {l^,S)i   namely, 


[S*;Z]    =  y   Z^  dS   =  0 
In  particular,   we  have? 


S* 


H+I   See  preface  for  references. 
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Corollaryf      There   exists   a  ■unique  harmonic   vector     Z* 
with  arbitrarily  specified      /^-periods   for  which 
Z^  =  0     on     S*. 

The  corollary  is   especially  aimple  in  not  requiring 
the  extended  definition  of  the   open  surface  periods;    all 
periods   are   functions  of  homology  class   alone   and  compati- 
bility is    automatic. 

As   an  example   of   the   above  corollary  consider  a 
vacuum  magnetic  field     B     in   a  domain  bounded  by  conductors. 
We  take     B     =  0      as   the   appropriate  boundary  condition* 
The  magnetic  vector*     B     will  be  uniquely  determined  by 
a  set   of  specified  fluxes     J,    =  ['7~,  fB] ,      or  by   the 
m.m.f.'s      [P^JB],      k  =   (l,2,..,,p)j      where    the     T~^^ 
constitute    a  basis   of   the   open  surface   group,    and   the 
I*,       are   a   reciprocal  basis   of    the   closed  ciirve  group. 
If  the      P,       boixnd  open   surfaces     ">",       in  the  exterior 
D     then  the  m.m.f.'s  may  be    interpreted   as   ciirrent 

A 

fluxes  given  by  uET",  ;J]   or  as  currents   ^jtl,   in  the 

/s 

external  circuits   f~\  • 

We  are  also  Interested  in  problems  for  which  the 
tangential  field  component   Z.   is  prescribed  on  the 
boundary.   Now  it  is  the  periods  on  a  representation 
^S;   or   [cj   of  the   C   group  which  must  be  specified. 

We  state  without  proofs 
Theorem  21   There  exists  a  unique  harmonic  vector  with 
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/s 


arbitrarily  assigned   C  -periods  (in  either  the  ^sj 
or   ^Cj   representation)  and  boundary  values  Z.   which 
are  local  surface  gradients  and  satisfy  the  compatibility 
condition  (14..2)  on  S  ,   but  are  otherwise  arbitrary. 
In  particular,  we  have, 

Corollary;   There  exists  a  unique  harmonic  vector  with 
arbitrary  assigned   (-^-periods  (In  either  representation) 
and  homogeneous  boundary  values  Z.  =  0  on  S  , 
Again  the  corollary  does  not  require  the  ejr tended  defini- 
tion of  the  open  periods. 

As  an  illustration  of  the  corollary  consider  a 
vacuum  electrostatic  field  E  in  a  domain  D  bounded 
by  conductors.   The  correct  boundary  condition  is  E.  =  0. 
The  electric  vector  is  determined  uniquely  either  by  a 
set  of  specified  potential  differences  ^v  "  'i'o  ~  t^,  jE], 

or  by  specified  external  charges  Q,  =  k;[S,  JE], 

Iff 

(k  =  l,2,.,.,q)|  where   the      C,       and     S,       denote   the 

reciprocal   bases   of   the   open  curve    and   closed   surface 

groups. 

Existence  theorems  for  the  inhomogeneous  differential 

systems,     ^ 

Turning  now  to  the  inhomogeneous  problems  (5»3)» 
(5.U)j  (5.1),  we  note  from  the  fact  that  the  solutions 
{S»5)t    (5*6)  exist  and  determine  the  general  solutions 
to  within  an  added  harmonic  vector,  that  a  meaningful 

-  60  - 


problem  may  be  expected  to  employ  the  same  sort  of 
boundary  condition  and  period  condition  as  for  harmonic 
vectors.   The  only  necessary  caution  is  that  it  now 
becomes  necessary  to  extend  the  definition  of  periods 
even  on  closed  cijrves  and  surfaces  to  make  allowance 
for  the  given  values  of   a'   and  ^  ,   The  values  of 
periods  must  be  assigned  to  fixed  representatives  of 
homology  classes  rather  than  to  the  classes  themselves. 

In  the  light  of  the  above  observations,  we  have: 
Theorem  3<   The  inhomogeneous  systems  (5«3)»  {S»k)f 
(5.1)  possess  unique  solutions  under  the  conditions  of 
Theorem  1  or  Theorem  2  with  the  compatibility  conditions 
referred  to  the  prescribed  fields  cr       and  ^  ♦ 

As  a  consequence  of  the  above  theorems,  we  shall 
prove  a  further  theorem  which  will  be  useful  in  the 
sequel. 
Theorem  ki      There  exists  a  unique  solution  of  the  system 

r  curl  curl  A  =  —  (div  ~^   =0) 

(5.8)     j  —  — 

(^  div  A  =  (T 

with  arbitrary  boundary  values     A^^      and  with  either 

arbitrarily  prescribed  periods      [C;A]      or   prescribed 

periods      [S;A3      compatible  with     cr  . 

The    accounting  of  boundary  data   to    see   whether   they 

seem   to  be   sufficient   to    specify   a   solution  is   fairly 

subtle.      The    equation      curl   curl   A   =  ^      is    akin   to    a 
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potential  equation  for  each  of  the  three  components  of 
A   (it  is  exactly  this  if  "^  =  0     and  <r  =0). 
However,  the  number  of  scalar  boundary  conditions  is  two 
(A^  is  a  two-component  vector) J  this  loss  of  one  boundary 
condition  is  necessary  in  order  to  satisfy  the  other 
differential  equation,   div  A  =  cT",   In  the  simpler  case 
where  only  one  "curl"  is  present,   curl  A  =  ^  ,  we  may 
impose  only  a  single  scalar  boundary  condition, 
A^  =  a  surface  gradient j   the  additional  "curl"  allows 
an  extra  boundary  condition. 

To  demonstrate  theorem  I4.,   we  set  curl  A  =  Z  and 
solve  in  turn 


(5.9) 

and 

(5.10) 


curl  Z  =  _^  ,       ( div  -  =  0) 

div  Z  =  0 


(div  Z  =  0) 


By  differentiating  Ax.  on  the  boundary,  (or  by  computing 
lim  -^    ^  A.    *  dx)     we  compute  the  boundary  values  of   Z  , 
and  obtain  the  periods   [2IJ2]   by  integration.   Since  A. 
is  given  as  single -valued,  these  values  are  compatible.   A 
unique  solution  of  (5«9)  therefore  exists  (compatibility 
with  respect  to  j^,  ^^  been  taken  care  of).   Now  we 
ttorn  to  (5.10).   The  boundary  values  Ao.   are  compatible 
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with  the  Irihoinogeneous  term  Z  by  ovr   construction  of 

Z  and  the  periods  are  assigned  either  In  the  form   [C|A] 

without  restriction  or  in  the  form  [SjA]  compatible 

with  gT"  by  hypothesis.  Hence  a  \anlque  solution  of 

(5»10)  exists,  and  the  theorem  is  proved. 

Theorem  I4.  has  application  to  the  system 

• 
curl  curl  E  =  -^i,J 

(^.11)     ^ 

div  E  =  q/K 

obtained  by  eliminating  B  from  the  pre-Maxwell  equations. 
Here   q  and   J  are  assumed  to  be  given  as  functions  of 
space  and  time.   The  boundary  value  E^  can  be  specified 
arbitrarily  and  also  the  total  charge  in  the  exterior 
domain  of  each  boundary  surface  (with  one  exception)  and 
this  will  determine  the  electric  field  uniquely.   Roughly, 
a  general  field  E   satisfying  (1.1)  is ' distinguished 
from  an  electrostatic  field  E   (q  =  0,  ()B/(^t  =0)   in 
that  E.   is  arbitrary  rather  than  merely  an  arbitrary 
surface  gradient.   The  specification  of  E  fixes  (}B/3t, 
but  not  quite  B|   some  initial  values  are  required 
for  this,  and  the  matter  will  be  discussed  in  Section  10. 
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6,   Orthogonal  Decompositions. 
In  working  with  the  system  (5.1),  It  Is  natural  to 
utilize  the  customary  viewpoint  of  linear  analysis  and 
to  write  a  solution  in  terms  of  the  decomposition  (5«7) 
into  harmonic,  Irrotational,  and  solenoldal  vectors. 
It  is  possible,  moreover,  to  make  refinements  in  the 
decomposition  so  as  to  mal:e  the  components  orthogonal 
in  the  sense  of  Hilbert  space  projections. 
Inner  product  formulas*   (Hodges'  Theorem.) 

We  Introduce  on  the  class  of  continuously  dlffer- 

l5 
entiable  vectors,  the  inner  product 

(6.1)  (X,Y)  =/  X.  Y  dV 

D 

A  simple  expression  for  the  inner  product  in  terras 

of  periods  and  boundary  values  alone  may  be  obtained 
when  one  vector  is  irrotational  and  the  other  solenoldal, 

f curl  X  =  0 

(6.2)  ^ 

I  dlv  Y  =  0 

There  exists  a  local  potential  such  that  X  =  grad  0, 
The  local  potential  may  be  taken  to  be  single  valued  by 
introducing  suitable  cross  cuts  in  D.   To  describe  the 
cuts  let  us  take  reciprocal  bases   P^,.,.,P'   in  ")  I^J 
and  XIi»***»ZIp   i^  ;^J»   (i«e.,   P^-   intersects 
J" ,     once  and  V^*      k  /^  j,   not  at  all).   Evidently 
the  potential  0     can  be  made  single -valued  when  we 
take  the   7~i   as  cuts  in  the  domain  D. 

15.   Note  that  we  are  not  dealing  with  a  complete 
Hilbert  space. 
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From  the  identity 
(6.3)        div  {0Y)   =   grad  0   'Y  +  j6   div  Y 
we  have  X  •  Y  =  div  {0Y)    -  f6   div  Y  =  div  (^Y)   and  upon 
integration  over  the  cut  domain 

f  X'Y  dV  =  f    j6y  'dS   +^      J      [jZJlY  .  dS 

where      [0]      Is   the    Jump  In     0     across     V",  «      Since 
[jZJ]    =   ^     X  •  dx      and     curl  X  =  0,      it   follows   that 

[^1      is  constant  on     5~i  »      if^l    =   CPt,JX].      Consequently, 

(6.1;)        (X,Y)   =  Xf  jZiY  .  dS  +  5^   [  Ti^jX]    iJl-^iY]      . 

The  period  notation  In    (6«[|.)   is  used  without   any  impli- 
cation that   the   values    are   constant   on  homology  classes. 
The    values   are   given  for   specific   curves       r~l       and   surfaces 
y ,  «      In  particular,    not   only   the   periods  but   also    the 
surface    integral  may   depend  on   the   choice   of  cuts.      The 
Integrals   on   the  right   depend  upon   the   particular   choice 
of   the   cuts     y~,  .      but   if  we   assume    the  boundary  condition 
Y     =  0  on  S    ,   we  have    (for   irrotational  X  and   solenoldal  Y), 


n 


(6.5)        (X,Y)   =  2^  [PkJX]    [ZkJY],  (Y^  =  0  on  S*) 

k=l 

If  on  the   other  hand     X^   =  0     on     S    ,      then  from 
Table   la,    ill,   we  have     ^   =  JZ$,  ,      constant,    on     S,  . 
Prom     /z5q  Y1   fS^*"^^    ~  °»      ^^   follows   that 
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Since     0,     -  0^   =  J„     X  •  dx,     we  have,   for     X     and     Y 

k  0         C^ 

satisfying   (6.2), 

(6.6)      (X,Y)   =  'fz  tCj^JX]    [S*}y]    ,  (X^  =  0  on  S*)   . 

The   results    (6.5)   and   (6«6)   together  constitute   a 
special  case   of  Hodge's   theorem  for  the   domain     D. 

In  summary,    for     X     irrotational   and     Y     solenoidal 
the    inner  product      (X,Y)      depends  on  boundary  values 
and  periods   alone    (6.i|.).      If,    in  addition,    either     X     or 
Y     satisfies    an  appropriate  homogeneous  boiondary  condition, 
the   inner  product   is   expressible    as   a  bilinear  form   in 
the   periods    alone,    for   one   case,    in   terms   of       C-periods 
(6.5)»   for   the   other,    in   terms   of     /?f-periods    (6.6). 
The   inductance   and  capacitance  matrices. 

If     Z      is  harmonic      Z   =  X  =  Y     in   (6.2)   and   satisfies 

either  boundary   condition     Z      ~  0     or     Z^   =  0,      then  the 

equations    (6.5)    and   (6.6)    express   the   Dirichlet   Integral 

/  Z      dV     in  terms   of  periods   alone.      If   the  boundary 

condition     Z     =  0      is   prescribed,    then     J   Z     dV     is 

^  D 

expressed  as   a  bilinear  form   in   the   periods   of  the   two 

representations  of  the  Jfl      group  J  hence,  by  taking 

Z  =  B,   the  magnetic  energy  ^  f  B   dV  may  be  expressed 

as  a  bilinear  form  in  cixrrents  and  fluxes.   Similarly, 

under  the  boundary  condition  2.  =  0,   the  Dirichlet 

integral  is  expressed  as  a  bilinear  form  in  the  periods 

of  the  two  representations  of  the   O  group J  taking 

Z  =  E,  we  have  the  electric  energy  ^jfj  ET  dV  expressed 

D 
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in  terms  of   charges   and  potential  differences.      In  each 
of  these  cases   the  periods  of  one  representation  of  the 
abstract   group  are   related  by  a  linear  transformation  to 
the  periods  of   the  other  representation*      In  the  first 
example   the   transformation  is    given  by   the    inductance 
matrix  and  in  the   second,   by  the   capacitance  matrix.      In 
either  case,   by  employing  the  matrix  it   is  possible   to 
express    the   Dirichlet   integral    as   a  quadratic   form  in  the 
periods  of   a   single   representation  of   the   abstract   group. 

Specifically,  we   consider  the   case  of   a  vacuum 
magnetic   field  in  a  domain     D     bounded  by  perfectly 
conducting  surfaces  upon  with     B     =  0.      Since   the 
prescription  of   one  class  of      ??t-periods   together  with 


the   condition     B      =0      determines     B     tiniquely,    the 
other  class  of    JT^  -  periods   is    also   determined.      The 
closed  curve   periods      [  P*,  |B]      are   determined  by  the 
c\ir rents      I,       in  the   external  circuits      H'v*     where 
[P,  ;B]    =  iJL   I,  .      The   associated  fluxes   are     J^  ~   [^^.jB]. 
To   exhibit   the    relation  between   the      I,       and    the     X  » 
we  represent     B      in  two  ways   as   the   sum  of   special 
solutions.      Under  the   condition     B      =  0     let     B.     be 
the  unique   solution  of   the   problem 

(6.7a)  f  ^k'^J^    "^  1^  \j  (J   =  !»•••»?) 

16.     The   Kronecker   symbol      b^.      is    defined  by 
_    To    ,   k  /^   j 
^kj  '   1^1    ,   k  =  J 
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t 
and     B  .      the  \inique    solution  for  which 

(6.7b)  CUkJBj]   =  \j      • 

The  field  B.   is  then  the  field  obtained  when  \init 
current  is  applied  in  the  one  external  circuit   P",  • 
Similarly,  B.  is  the  field  for  which  unit  flux  exists 
across  ]^.  but  the  flux  is  zero  across  every  other 
independent  open  surface •  For  an  arbitrary  solution  B 
specified  either  by  its  currents  I^,Ip,.»«,I   or  by 
its  fluxes  5-|»5p>»«»»5  f      we  may  then  put 

B  =  I^B^  +I2B2  +  ...  fXpBp 

(6.8)  ^   ,   ^   ,        ^   , 

This   is   a   specific   exemplification  of  the    statement 
that    the   linear   space   of  harmonic   vectors   in     D     is 
finite   dimensional.      In  particular,  for  each     B.     we  may 
calculate    the   fluxes   and  obtain 

(6.9)  B^.   =  5;^  L^^  B^   ,  (J=l,...,p), 
where     L.,       is   defined  by 

(6.9a)  L^^  =  tUk'^j^      ♦ 

Hodge's   Theorem   (5*6)   gives 
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where  by  taking  B,   and  B,   in  turn  as  X  we  see  that 

the  matrix  L.,   Is  symmetric.  Prom  (6.8)  and  (6«9)  we 

P.  .     p  , 

have     B  =  >      ^ ^-i^iir^v     ^^^  obtain  the  fluxes  by 

j=l  k=l  ^  ^^  ^ 

(6.10)      ik-fll-jk^j   J 

J  •*• 

we  see  then  that  L.,   is  the  classical  Inductance  matrix. 

The  magnetic  energy  now  may  be  given  the  simple 
representation, 

D 


I  J  ^  dV  =  ^  1^  [  r,,Bl  [Z,<B1 


(6.11) 


=  I IZ  Ik  i„  =  E  Ik  ^"-  ^ 


from  which  it  follows  that  the  matrix  is  positive  definite. 
If  the  basis   P,   written  as  a  vector,  is  transformed 
to  a  new  basis   P*   by  a  matrix   (M),   then  YI  ^^  trans- 
formed  to  a  reciprocal  basis  Z-.     ^J  ^   related  transformation 

r'  =  (M)  r  ,    r'  =  (M*)-^r: 

where   (M  )    is  the  reciprocal  transpose  of  M, 
Accordingly  we  have 

i'  =  (M*)"^  i  =  (M*)-^  (L)  I  =  (1^-1(L)(M-^)  l' 
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so  that   (L)  transforms  according  to  the  rule 
(L)  >  (M^)"^  (L)  (M"^) 

The  induction  matrix  is  \iniquely  determined  by  the 
domain  D  and  by  the  choice  of  a  basis  for   vt ,   that 
is  a  choice  of  external  circuits  carrying  assigned 
cvirrents   I,   or  of  surfaces  with  assigned  fluxes  J,  • 
It  is  important  to  remember  that  to  express  the  energy 
in  the  simple  form  (6»ll)  it  is  necessary  to  take 
reciprocal  bases  in  defining  the  currents  and  fluxes. 

In  direct  analogy  with  the  preceding  case  there 
is  a  symmetric  positive-definite  capacitance  matrix 
relating  the  potentials  to  the  charges  if  the  Imrmonic 
vector  is  E   and  the  homogeneous  condition  E.  =  0  is 
imposed.   The  e,m»f. 's  and  charges  are  defined  on 
reciprocal  bases  by  ^v  "^  i^>  ~  )^0  ~  t^j^jE]   and 
\  ~     *^tSj^JE]  {  (k  =  l,2,,,.,q).   These  are  related  by 
the  classical  capacitance  matrix 

(6.12)       ^=^G^jVj  (i  =  l,2,...,q) 

where  C.  .   is  defined  in  the  manner  of  the  preceding 
example.  The  electrostatic  energy  is  given  by 
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iKJ^ 


(6.12a) 


'^  k=l     "^ 


^k  «k3  '3 


Since  we  asstirae  for  the  total  charge   of   space 
we  may  write    (6.12a)   in  the   ayrametric  form 


\ 


=  0 


(6.12b)  \\     E^  dV  =  I  7"  ^. 


Sc 


Orthogonal   decomposition  of  vectorsa 

Hodge's   theorem  yields   a  number  of  useful   refine- 
ments of   the   decomposition   (5«7)   of  an   arbitrary  vector 
into   irrotational   and   solenoldal   components©      First  we 
observe   that     X     and     Y     are   orthogonal,      (XjY)   =  0, 
if  either  of   the   conditions   in  the  following  table   is 
satisfied. 

Table   II 
Orthogonality  Conditions  for   Irrotational 
and   Solenoidal  Vectors ^      curl  X  =  0^,    div  Y  =  0. 


.-» 


X^  =  0  on  S  o   [C|X]  or  [S|Y]  =  0 


Y^  =  0  on  S^^o   C  P|X3  or  EHjY]  =  0 


Prom  Table  lib  it  follows  that  the  harmonic  vectors  L 
with  L  =  0   are  orthogonal  to  the  harmonic  vectors   M 


n 


with      [ PjM]    =  Oo      From  Table    lla  harmonic   vectors     N 
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with   N^  =  0  are  orthogonal  to  harmonic  vectors  M 
with   [SjM]  =  0.   The  three  classes  of  harmonic  vectors 
satisfying  the  conditions  of  the  following  table  are, 
then,  mutually  orthogonal. 

Table  III 
Orthogonality  Conditions  for  Harmonic  Vectors 


L 

M 

N 

L„  =  0 

[  PjM]    =   [SjM]    =  0 

Nt   =0 

The  value  of  Table  III  lies  in  the  fact  that  the 

decomposition  is  complete  as  well  as  orthogonal.   In 

other  words,  every  harmonic  vector  may  be  assigned  a 

unique  representation  as  a  sum  of  orthogonal  harmonic 

vectors,  one  from  each  class.   Let   H  be  an  arbitrary 

harmonic  vector  which  takes  on  the  boundary  values 

H  ,H,   and  has  the  periods   [PjH],  [SjH].   Prom  the 
n   u 

existence  theorems  for  harmonic  vectors  we  can 
determine,  in  turn,  the  three  unique  vectors 

La     \  =  0  ,     [  PjL]  =  [  Pjh] 

NJ      N^   =  0    ,       [SiN]    =   [SjH] 
MJ      M  =  H-L-N      . 

Plainly,      H=L+M+N     Is   the  required  decomposition. 
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In  the  same  way,  by  applying  the  existence  theorems 
for  harmonic  vectors  we  obtain  from  Table  II  an  orthogonal 
refinement  of  the  decomposition  (5»7)  of  an  arbitrary 
vector  Z  Into  Irrotational^  harmonic,  and  solenoldal 
components!  namely,  components  satisfying  the  conditions 
of  the  following  table. 

Table  IV 
Orthogonality  Conditions  for  Irrotatlonal« 
Harmonic  and  Solenoldal  Vectors, 


X  (=grad  f6) 

H 

Y  (=curl  A) 

curl  X  =  0 

div  H  «  0 

dlv  Y  =  a 

Xt  -0 

curl  H  =  0 

Y^  =  0 

[C}X]  =  0 

[J    JY]  =  0 

Again,  this  decomposition  is  not  only  orthogonal  but 
also  complete!  every  vector  Z     dan  be  uniquely  expressed 
as  the  sum  of  three  vectors,  one  from  each  column.   The 
headings  c\arl  A   and  grad  0     are  indicative  of  the 
fact  that  single-valued  potentials  exist  under  the 
conditions  stated.   The  middle  column  can,  of  courae, 
be  further  subdivided  as  in  Table  III.   This  is  the 
orthogonal  refinement  of  (5»7)  that  was  previously 
referred  to. 

To  prove  the  completeness  of  Table  IV,  we  start 
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with  the  non-orthogonal   decomposition   (5«7), 

Z  =  Hq  +  grad  0^  +  curl  Aq      • 

There   exist  unique  harmonic  vectors,     H,      with  the   same 
tangential  boundary  values   and     {cj      periods   as     grad  /ZJ^, 
and     E^     with  the  same  normal  boundary  values   and     )  5~  j 
periods   as      curl  Aq.      The  period   and  boundary  data  are 
clearly  compatable.      Setting     H   =  Hq   +  H,    +  H^    , 
grad  jZJ  =  grad  0^   "  %»      °^rl  A  =  curl  Aq   -  H^,      we 
©btain  the  decomposition  in  the   desired  form* 

It  will  be   convenient   also   t©  have   this   result  in 
the  form  of 

Theorem  5*      Any  arbitrary  vector     Z     can  be  written  in 
the  form 

(6.13)  Z  =  curl  A  +  H  +  grad  0 

where     f6  =  0      and     A.    =  0   on  the   boundary. 

Theorem  5  follows   from  the   fact   that   the   first 
column  in   IV  describes    the   gradient   of   a  potential 
function  which  is   constant  on     S        and  may,    therefore, 
be   taken  as   zero    there.      At   the   same   time,    the   third 
colvimn  in  IV  describes   the   curl   of  a  vector  potential 
A     where      A.      is   a   surface  gradient    and 

/      A4.»dx  =  0  (k=l,...,p)« 

P 
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Consequently,  by  subtracting  a  suitable  hannonic  vector 

t  t 

from  A,  we  can  obtain  an  A   for  which  A.  =  0  « 

We  have  used  the  existence  theorems  of  the  last 
section  to  prove  completeness  of  the  orthogonal  decomposi-f 
tions  of  Tables  III  and  IV.   A  more  common  procedure  is 
the  reverse  J  namely,  to  enlarge  the  space  of  differentiablo 
vectors  X  and  Y  to  obtain  a  complete  Hilbert  space, 
appropriately  generalizing  the  differential  relations 
(6*2)  and  the  bo\Andary  conditions  so  that  they  do  not 
imply  smoothness.  The  existence  and  uniqueness  of  the 
orthogonal  decompositions  then  follow  by  application  of 
the  projection  theorem  in  Hilbert  space.  A  proof  in 
this  way  of  the  existence  theorems  of  the  last  section 
is  trivial  in  the  class  of  Hilbert  space  elements,  but 

additional  effort  is  required  to  verify  actual  differ- 

17 
entiability  of  the  harmonic  vectors   • 

Further  inner  product  formiilas. 

Since  we  shall  be  dealing  with  inhomogeneous  systems, 

it  is  not  sufficient  to  work  with  the  representations  of 

inner  products  given  by  (6.5)  and  (6.6).   For  the  system 

inhomogeneous  in  Y 


curl  X  =  0 

(6.14) 

div  Y  =  (T 

17,   Cour ant -Hilbert,  Methoden  der  Mathematischen 

Physik,  Berlin,  Springer  (1937),  Vol.  II,  Chap.  VII, 
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equation  (6«3)  yields 


(6.15)     (x,Y)  =  -  J  jzJ  srdv  +  J    izJy  .  ds  +  XI  f  Pk^XlCZk^"^^ 


D  S' 

For  the    system   inhomogeneetis   in     X, 


(6«ad) 


curl  X  =  7 
div  Y  =  0     , 


(div    -    «  0) 


matters   are   somewhat  more  complicated* 

To  obtain  a  single-valued  vector  potential     A 
with     Y  =  curl  A,      cuts   are  made    along  the   ©pen  arcs 
C, ,«»»,C    •      To  perfoiTO   the  necessary  integrations,    the 
cuts     C,       are   enclosed   in  narrow  tubes     T,     (k  =  l,»».,q). 
Figure  10.      Upon  integrating  the   identity 


(6.17)  div    (A'^X)   =  X -curl  A   -  A  •  curl  X 


Figure  10. 


A  domain  cut  by  an  open  arc      C,      for  the 
purpose  of  obtaining  a  single-valued  vector 
potential.    " 
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we  obtain,    excepting  small  elements^ 


(X,Y)   =   {A,jz^)   +  J     A«X  •dS  +  211    f   AaX  -dS 

«  k=l 

First  we  note  that 

/a  X  X  •  dS  =  fk^  «  X^  •  dS   . 

If     T,       Is   permitted  to   shrink  t©      C^,      the   integral   of 
A     aroTind  the   tube   approaches  the   constant      (S,  JY]      (by- 
Stokes'    Theorem),    and  since     X     approaches   a  value 
depending  enly  on  the   limiting  position  on     C^ ,     we   find 

/   AxX-dS  >  [S^IY]    [Cj^jX] 

We  note   that      y  A  *  X  •  dS     converges    since   the   singularity 
in     A^      is   integrable.      Finally  we  have 


(6.18)      (X,Y)   =   (A,jr.)    +/    A.'^X.    .dS  +  2_   [^}Y][Cj^jX] 

S*  ^=1 


Each  of  the   terms   on  the  right    in    (6.l5)    and    (6.18) 
depends   on  the   choice   of   cuts.      The   field      A     may  be 
singular   at   a  cut   and     0     cannot   in  general  be    extended 
in  single-valued  fashion  across  cuts. 

If     €r=  0      and     ^  =  0,      (6.l5)   and   (6.18)   combined 
yield  a  relation  between  siirface   fields    and   periods, 
namely. 
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7rc-    "ills   we    0  ai:   oz'air:   a.  usef-il   iienti-y   i:iTolTlr^ 
rwc    arbl-rsLr-   rjLT-acs    gridift--s,    say,     l_      ar.!     X_ . 

(6.19)    "  A_^  x_  -is  =  ^^__  cTvixitrLsA:  -  cr^jx] :"".,}!]  . 


/N 


,* 


Herfl       r,       is   -i^e    clased  0":irTe   ©r.     S        viich.  bo-JirLds     >  .. < 

£  — - 

fr'f  g    forrr^a   1=    a   or"  serience   of   "is   ob  ssrraTiozi   "liia.":    an 
BTz LzT s.r-r   s-JLrfac*    sradien"      X     car.   ce   represented  &a   the 
dirferencs    of    TaO    rectors,      X        and     X   ,      wnere 
z-jLTl   X      =   :      in      3      and      c-.Lrl   X     =  C      in  the   ccnplenent 
of      Z        in  this   vay   an     X_      can  re    oonstrticted  with   eJll 
its    surface    z-ericis    ircitrary   and    different   fr«n  zer«)« 
I:   "ce    s--:re,    C6.19)    can  ce   obtained  siore    directly  by 
integrating  the    .iacal)   identity 


c-^1    j;  T3V   -   ^  Tg^  ;    =  2  ^3^  «  73! 


over 


surface  with      P^,  1"^,    ''^  =  1,..,,?)      introduced   as    cuts, 

Fortrnla      6.19)    i3    a   tvo-di^ensienal    's-^face) 
analogue   of   the   rrevirus    three-dinensienal   inner   troduct 
ferrrilAS.      Ihe    3-:Lrf&ce    gradient,      A.,      has    its    surface 
curl   zero    and   the    oeabination     X^^  *  ^     vhene      X_      is    a 
surface    ss'a.dient    is    easilv    seen  to   ha.Te    zero    surface 
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divergence,   properly  defined.      The    analogy  is    somewhat 
hidden  since    (6.19)   applies   to    a  closed  2-inanifold, 
whereas   the   three-dimensional  formulas  referred   t© 
open  3-manifold3  with  boundary.      The   analogous  result 
w©\ild  be   like  formula   (9.3)  with  the   summations   extending 
over  the  full  groups,    e.g.,   by  considering  the  ©pen 
periods   as  closed  by  some   sort  ®f  boundary  idontiflcati©n« 
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7.   Two-dimensional  Harmonic  Fields. 

18 
Cylindrical  Problems 

The  theory  of  two  dimensional  harmonic  vectors  Is 
particularly  elegant.   The  loss  of  one  dimension  leads 
to  a  homology  theory  which  Is  not  merely  a  specialization 
of  the  three-dimensional  homology  but  has  its  own 
intrinsic  symmetries »   For  example,  the  homology  properties 
of  closed  surfaces  and  closed  curves  are  unrelated  in 
three  dimensions  but  the  closed  curves  in  the  two-dimensional 
case  are  the  counterparts  of  both  the  closed  curves  and 
the  (projection  of)  closed  surfaces  in  the  three-dimensional 
case.   This  fact  is  reflected  in  duality  relations  which 
virtually  halve  the  work  of  formulating  the  theory  of 
plane  harmonic  fields  and  their  kin. 
The  planar  theory.   Duality  in  the  plane. 

As  before,  our  attention  is  devoted  to  the  solution 
of  problems  of  the  form 

div  Z  =  0" 

(7.1)         j 

curl  Z  =  T 

where  the  divergence  and  curl  are  now  operators  in  the 

plane,  and  the  scalar  fields,  c,  t  are  assumed  to  be 

19 
given  data  of  the  problem   .   In  the  plane,  however. 


18.  This  section  may  be  omitted  or  postponed  without 
disturbing  the  continuity  of  exposition. 

19.  Since  the  planar  curl  has  only  one  component  it 
may  be  treated  as  a  scalar. 
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there  exists  a  duality  between  the  divergence  and  cupI 
which  permits  a  considerable  simplification  of  the  problem. 
To  each  plane  vector  field  Z  we  may  associate  a  dual 
orthogonal  field  Z,  where,  if  Z     is  given  component- 
wise in  a  cartesian  representation  by  Z  =  {Z,,Zp)   then 
Z     is  defined  by 

A  A 

(7«2)       Z^  —  Zq  j    Zo  ~  •"  Z-|   • 

With  this  convention  the  duality  of  the  vector  operaters 
is  exhibited  by 

curl  Z  =  div  Z 
(7.3)     < 

div  Z  =  -  curl  Z  • 

The  divergence  theorem  in  the  plane  for  Z  will  then 
prove  to  be  the  same  as  Stokes'  theorem  for  the  dual 

A 

vector  Z. 

We   shall   let     ^     denote   a  plane   domain   since 
this  notation  will  be   consistent  with  the  notation  for 
^pen  surfaces   in   three-dimensional  problems »      For  a 
simply  connected  plane   domain     ^     with  closed  bound- 
ing curve      r*,      the   integral   theorems   take   the   form 

f  div  Y  dS  =  ^  Y  •  n  ds 

n  r 

y*  curl  X  dS  =  ^  X  •  t   ds 

Z  P 
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(T.ll.) 


/^ 


where  n  =  t  denotes  the  exterior  normal  to   P.   The 
extension  of  the  Integral  theorems  to  more  complicated 
topologies  is  greatly  simplified  by  the  duality  relations, 
To  this  end  let  us  investigate  the  homology  properties 
of  plane  domains,  ^  ,/ 


Figiire  11.   Two-dimensional  homology. 

Let  ^  be  a  plane  domain  bounded  by  a  finite 

number  of  simple  closed  c\irves   fo>r^»»»«>Pp  oriented 

so  that   P'q  +  Pf  +  ...  +  n*  '-^  0.   The  closed  curves 

r^3»"  »\~^     serve  as  a  basis  for  the  closed  curve 
1'      r 

homology  group  ]  P/  in  Xl»   ^°^  ^  basis  of  the  open 
arc  homology  we  may  choose  arcs  Q.fC^f**^^     where 
each  arc  G,   connects  a  point  Pq  on   Pq   to  a  point 

The  analogy  between  the  plane  homologies  and  the 
three-dimensional  homologies  is  easily  visualized  in  a 
cylindrical  domain  D  which  has  ^  as  its  plane 
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cross-section.   The  closed  curves  in  ^  are  the  pro- 
jections of  both  the  boimded  closed  curves  in  D  and 
of  the  closed  surfaces  which  have  infinite  cylinders 
as  a  basis.   However,  the  closed  curve  group  in  D 
has  an  additional  basis  element  which  may  be  taken  as 
an  infinite  line  perpendicular  to  ^«   In  ]^,   the 
dual  group  to  the  closed  curves  are  the  open  arcs,  and 
these  correspond  to  projections  of  both  the  open  arcs 
and  the  open  surfaces  which  have  as  a  basis  infinite 
cylindrical  ribbons.   In  addition,  the  open  siorface 
group  in  D  has  ^  itself  as  a  basis  element. 

The  results  relating  periods  and  boundary  values 
for  an  Irrotatlonal  X  are  formally  identical  with  the 
results  in  three  dimensions  (Table  la.  Section  3»  page  hh) o 
Moreover,  it  is  completely  unnecessary  to  put  the  results 
for  a  solenoldal  vector  Y  (Table  lb)  in  planar  form. 
It  is  sufficient  to  write  the  table  for  an  Irrotatlonal 
vector  ±h   terms  of  the  dual  of  Y  through   dlv  Y  =  -  curl  Y, 

Y  =  -Y.*      We  have  only  to  observe  for  a  given  curve  with 

/\ 

unit   tangent      t      and  normal     n  =  t      that 

^t   "  ^n   *  ^n  "   "  ^t 

to  obtain  relations  between  the  periods  Involving  the 
boundary  normal  component  of  a  vector  and  those  involving 
the  tangential  normal  component  of  the  dual.   In  particular, 
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the  ^-periods  of  Z  are  given  by 

'  k 

J     Z     ds   =  -  J     Z.  ds  =  -  [C,  iZ]   , 

and  the   c  -periods  by 

/  2^  ds  =  -  >"  Z,  ds  =  -  [  r^jz] 

/  Z^  ds  =  [C^iZ]      . 

The  existence  theorems  also  dualize  quite  nicely. 
As  before,  the  solution  of  the  system  (7«1)  is  decomposed 
into  solutions  of  the  systems 


(7.5) 


(7.6) 


(7.7) 


C   div  Z  =  0 
(  curl  Z  =  0 

div  Z  =  cr 
curl  Z  =  0 

div  Z  =  0 

curl  Z  =  r 


For  (7.6)  we  have  the  special  solution 
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z  =    v/zi 

(7.8) 


f^  =  -^  J   <^^og  r     dS  , 


2it 

z 

and  by  putting   (7»7)   into   the   dual   form  we  have   the 
solution  of    (7.7)   directly  from   (7.8)   in  the  form 

2  =    V  \lf 

(7.9)  [ 

r 

The   functions     f6     and     \1;     correspond   to   the   concepts   of 
velocity  potential   and   stream  function  in  steady 
incompressible  plane   flow#      In   analogy  to    (5.7)  we   see 
that   an  arbitrary  plane  vector  may  be  written  in   the  form 

(7.10)  Z  =  H  +     ^^  +      V\l/      , 

as   the    sum  of   a  harmonic   vector,    a  gradient    and   the    dual 
of   a  gradient. 

By   virtue   of   the   representation    (7«10)   the   diffi- 
culties   in  formulating    a  well-posed  problem   for   (7.1) 
are    relegated  to    the  harmonic   parto      Let  us    then   consider 
harmonic   vectors   in   the   domain     ^«      In  the    simplest   case 
where     ^     is   the   Interior  of   a  simple   closed  curve      P^^ 
there    is   a    single-valued  potential     0     with      Z  =     V0 
and      S7   0  -  0,      The  Laplace   equation  for     fi     has   a  unique 
solution  if,    for  example,    either     (6     is   given  on      f"'^ 

-  85  - 


or  d^/^n     subject  to  ,P^   4^  ds  =  0.   It  follows  that 
a  ixnlque  Z  exists  either  with  prescribed  boundary 
values  Z^  subject  to   ^\^  Z^  da  =  0,  or  boundary 
values  Z^  subject  to    ^^  Z^  ds  =  0.   It  is  no 
accident  that  the  compatibility  relations  for  the  two 
problems  are  similar.  The  problem  for  harmonic  Z 
with  tangential  boundary  values  is  identical  with  the 

A 

problem  for     Z     with  normal  boundary  values.     More 
generally  when  the   boiindary   is  given  by 

a  haifmonic   vector  will  possess   two  kinds  of   closed 
curve   periods; 

[r;Z]    =o/  Z^  ds    (=  ^Z^  ds) 

(7.11) 

[  PlZ]    =  ^  Z^  ds  =  -  ^  Z^  ds      . 

Tangential  boundary  values      Zx.     may  be   freely  prescribed 
on      n^     subject   to   the   compatibility  requirement 

(7.12)  [P^IZ]    =   JZ^    Z.    ds   =  0      J 

similarly,    the  prescription  of  normal  boundary  values 
is  subject   to   the    compatibility  condition 

(7.13)  [r*jZ]   =  -   /    Z_    ds  =  0     . 

r 
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Homogeneous  boundary  conditions  yield  open  arc 
periods  defined  as  functions  on  the  open  arc  homology 
group,  namely, 

[CJZ]  =  /z^  ds  (Z^  =  0  on  P*  ) 

C 

(7.1U) 

[CJ2]  =  -  /  2^  ds  (Z^  =  0  on  P*) 

C 

There  is  no  longer  any  distinction  between  the 
C/    group  and  the  nl    group {  yet  a  shadow  of  the 
distinction  remains  in  the  designation  of  the  -^-periods 
of  Z  by  [  PjZ]   and  [0|2]   and  the   (^f^-periods  by 
[  PjZ]   and   CC|Z],   but  the  ^-periods  of  Z  are 
C/-periods  of  Z  and  vice  versa. 

Quite  as  before,  if  inhomogeneous  boujidary 
conditions  are  given  the  concept  of  open  period  must 
be  extended  and  the  values  of  the  periods  on  a  homology 
class  are  defined  in  terms  of  the  value  on  a  fixed 
representative  and  the  known  boundary  values.   With  this 
extension  of  the  concept  of  period  we  have  the  plane 

analogue  to  Theorem  Ij 

t 
Theorem  1  ?   There  exists  a  unique  harmonic  vector  Z 

in  ^  with  arbitrarily  prescribed  AZ-periods  (either 

[PjZ]   or   [C}Z])   and  boundary  values  Z   on   P 

subject  only  to  the  compatibility  condition  (7.13)o 

The  dual  of  this  theorem  is  the  analogue  of  Theorem  2, 

r 
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namely J 

Theorem  2    i     There   exists   a  unique  harmonic      Z     in     > 
with   arbitrary      O-periods      ([PjZ]      or      [CjZ]  )      and 
boundary  values     Z^     on      P        subject   only  to    (7»12)» 
Again,    by  extending  the  concept  of  period  to   the 
inhomogeneous    equations    (7«1)   so   that   the    values  of 
the  periods  on  a  homology  class  are   calculated  from 
the   value   on  a  fixed  representative,    the  known  fvmctions 
0^      and    TT,      and  the   prescribed  boundary  values  we 

obtain  the   analogue    to   Theorem  3  J 

t 

Theorem   3    i      Th©    system   (7»1)   possesses   a  unique 

solution  with  arbitrarily  prescribed     //l-ipevloda 
and  normal  boundary  values     Z        stibject   to 
,>£_^  Z     ds   =  */    0'dS     or  with   arbitrarily  prescribed 
^-periods    and  tangential  boundary  values     Z.      subject 
to     cX^     Z.    ds   =  y  r  dS   . 

r*    t       ^ 

Hodge's   theorem  for  the   inner  product   of  tw© 
vectors,    an   Irrotatlonal      X     and    a   solenoidal     Y, 


(7.15) 


follows   as  before   with   the   arcs     C^,C^,. ,• ,0^     taken 
as  cross-cuts   in     ^,      For  homogeneous   normal  boundary 
conditions  we  obtain 

(7.16)    (X,y)   =  /  X'Y  dS  =  "TZ  [r*jX]    [Cj^jY]      , 

(Y„  =  0   on    r*)   • 
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as  the   analog  to    (6,5).     The   analogue  to    (6,6)  is  simply 
the  dual  to    (7 •16)   and  Is   obtained  directly  from  the 

A      A 

observation  that      (X,Y)  =  (X,Y)     in  the   form 


(7.17)   (X,Y)  =  H:  [C.  JX]    [r*>Y3    ,  (X^  =  0  on   P*). 


k=l 

This  duality  has  a  curious  consequence  Tor  harmonic 

electric  and  magnetic  fields.   In  the  planar  case,  the 

inductance  and  capacitance  matrices  are  essentially 

inverse,  in  contrast  to  the  three-dimensional  case  in 

which  the  two  are  completely  unrelated. 

The  inductance  matrix  for  harmonic  B   in  J" 

under  the  condition  B   =  0   is  defined  by 

n  '' 

where  |x  I^^  =  [  Pj^jB]  ,  J  =  [Cj^jB],  Similarly,  the 
capacitance  matrix  for  harmonic  E  with  E.  =  0  is 
defined  by 

r 


where  \  =  "    K  [  P*;E]   and  V^  =  [Cj^jE].   Since  the 
problem  with  homogeneous  tangential  conditions  on  E 
is  equivalent  to  the  problem  with  homogeneous  normal 

A 

conditions  on  E  we  make  the  identifications 
\  =  "'^f^'k  ^^^  =  -  ^^   ^k 
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from  which  It  follows  that 

(7.18)  (Cj^)  =  o^(Lj^)-l 

where  c^  =  I/KTm-.   In  the  three-dimensional  case  the 
(^  and/'Lgroups  generally  have  different  dimensions 
so  that  even  the  orders  of  the  two  matrices  differj 
in  the  planar  case  the  inductance  and  capacitance 
matrices  are  practically  inverses. 
The  cylindrical  theory. 

The  term  two-dimensional  may  be  considered  more 
generally  to  apply  to  cylindrical  problems  rather  than 
purely  planar  problems.   In  other  words ,  a  two-dimensional 
problem  is  one  in  which  the  fields  and  the  attendant 
conditions  are  independent  of  one  cartesian  coordinate, 
but  vectors  are  permitted  a  non-vanishing  component  in 
the  third  direction.   Let  x  and  y  denote  the  relevant 
cartesian  coordinates  and  let  z  be  the  ignerable 
coordinate.   In  problem  (5'l)  we  set 

z  =  z'  +  z*,   ~r-  =~r-«+-r^ 


where  Z  ,_-_       are  the  components  in  the  direction  of 
the  z-axis.   The  problem  is  easily  reduced  to  the 
solution  of  the  problems 

f  div  Z  '  =  cr 

(  curl  Z  '  =  ^^ 
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and 

(7.'20)       grad  Z*  =  -  H'  (°^rl  Z!'=  ^^  • 

The  cylindrical  problem  may  then  be  decomposed  into  the 
plane  problem  (7«19)  and  the  elementary  problem  (7<»20)« 
The  cylindrical  homology  does  differ  slightly  from  the 
plane  homology.   If  we  let  ^Z  represent  the  plane  cross- 
section  of  the  cylinder  it  is  clear  that  the  cylindrical 

^-group  is  Isomorphic  to  the  two  representations  of  the 
plane  group.   On  the  other  hand  the  cylindrical  ^T^group 
has  an  additional  basis  element  which  may  be  taken  as 
y     itself  in  the  open  surface  representation  or  as  a 
single  infinite  line  perpendicular  to  ^  in  the  closed 
curve  representation.   Coupled  with  this  additional  basis 
element  we  have  a  trivial  period  condition  for  problem 
(7»20)  which  requires  the  prescription  of  either   «/  Z  dS 
or  Z  (Xq)   for  arbitrary  Xq   in  Yl»      However,  for  the 
single-valuedness  of   Z   in  (7«20),  the  given  function 

— ^   must  satisfy  beyond  curl  j^.  ~  ^y      '^®  stronger 
condition 

(7.21)  /  ~'  da  =  0  . 

r— I  — n 

for  all   closed  ciirves      f^     in     ^   . 

The   cylindrical   electromagnetic   equations* 

The   cylindrical   problems  for  the   pre-Maxwell   system 

(1.1)  may  easily  be   decomposed   according  to   the    above 
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scheme.  For  example,  the  pre -Maxwell  equations  may  be 
solved  first  for  B  and  then  for  E  by  employing  the 
eq\xations 


curl  B  =  tiJ  ,      grad  E  =  4| 


'v « 


div  b'    =  0 


(7.22) 


grad  B     =  -  |xJ  curl  E     +  4|-  =  0 

(^J^  ds   =  0)  div  e'    =  q/^: 


or  first   for     E     and  then  for     B  by  means   of   the 
equations 

V  E,    -  yi  -J^     ,  ^  =  graa  E 

grad   curl  E     =  iJ,  ^^  ,        H^  =  -  curl  E 


div  E     =  q/i^ 


With  this  decomposition,  the  existence  theorems  developed 
in  Section  10  may  be  immediately  interpreted  in  the 
cylindrical  case.   These  methods  have  immediate  appli- 
cation to  the  cylindrical  pinch  problem  treated  in 
Section  10. 
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8.  Variational  Methods   . 

We  shall  now  show  how  to  construct  a  series  of 

variational  problems  with  unique  solutions  which  are 

exactly  the  solutions  of  the  differential  problems  of 

Theorems  l-3<>   This  will  be  a  convenient  way  of  looking 

at  these  questions  In  a  certain  class  of  fluid  magnetic 

stability  problems  which  we  shall  treat  in  a  later 

report  in  this  series.   In  each  of  the  existence  theorems 

a  vector  is  uniquely  determined  by  differential  relations, 

dlv  Z  =  cT"  and  curl  Z  =  j^.*   ®^^  ^7   given  boimdary  values, 

Z^     or  Z  ,   and  certain  period  specificationsa  By 

dividing  the  data  which  xmiquely  determine  Z  Int©  two 

categories  to  be  labeled  "ac3missible"  and  "variational" 

(their  precise  specification  will  be  given  later),  we 

shall  characterize  Z  uniquely  as  that  vector  of  the 

class  having  the  same  admissible  data  as   2  which 

minimizes  a  particular  functionalo   In  other  words ,  the 

vector  Z  will  appear  as  the  unique  minimizing  vector 

for  a  variational  problem  in  which  the  variational 

conditions  play  the  role  of  Euler  equations  and  natural 

boxondary  conditions »   A  simple  example  of  a  problem  of 

this  type  is 

2ai;      Cf.   Courant-HllbertI      Methods  of  Mathematical 
Physics,    Interscience    (1953),    Chap.    IV. 
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admissible  I   curl  X  =  0  ,  X^  given 
variationals   div  Y  =  0   ,   [SjY]  given  • 

The  manner  in  which  the  functional  appropriate  to  this 
problem  is  obtained,  will  be  described  below. 
Formulation  of  the  variational  problem* 

The  data  for  Z  are  divided  by  the  following  table 
into  two  natural  categories  according  to  their  relevance 
in  the  period  conditions  (cf.  Table  I),  regardless  of 
their  eventual  designation  as  admissible  or  variational! 

Table  V 
Division  of  Boxindajry  and  Period  Data 
into  Curl  and  Div  Classes 


a. 
b. 


curl   Z,  Z.,  [n;Z],  [CJZ] 


div  Z    ,  Z^,  [HiZ],  [S;Z] 


For   the   present,    the   table    is    to   be   considered  merely 
as   a  mnemonic   device    in  which   tangential   boundary  values 
and   curvilinear  periods   are    associated  with     curl   Z, 
while  normal  boundary  values   and   surface   periods    are 
grouped  with     div  Z.      Corresponding   to   this   division 
of  data  the   data  in  a  well-posed  problem  for     Z      are 
assigned  to    tw®   classes   of  vectors,    the  curl   class 
consisting  of  vectors     X     which  agree  with     Z     in  those 
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data  which  fall  in  line  a  of  the  table,  and  the 
div  class  of  vectors  Y  which  agree  with  Z  in  that 
part  of  the  data  falling  in  line  b.   In  the  example 
above,  the  curl  class  is  taken  as  the  admissible  class 
and  the  div  class  as  the  variational  class » 

In  general,  corresponding  to  every  set  of  data 
which  determines  a  unique  vector  according  to 
Theorems  1-3,  we  shall  have  two  distinct  variational 
problems  depending  upon  which  one  of  the  curl  class  or 
div  class  is  called  admissible  and  the  other  variational. 
The  two  related  problems  are  called  conjugate.   Further, 
there  is  a  precise  duality  between  the  categories  (a) 
and  (b)  of  Table  Vj   in  all  the  existence  theorems  and 
orthogonality  relations  a  statement  remains  true  if  curl 
is  replaced  by  div,  tangential  condition  by  normal,  etc. 

There  remains  the  problem  of  assigning  a  variational 
fionctionj  that  is,  a  functional  which  is  to  be  minimized 
on  the  admissible  class.   Given  an  association  of  the 
data  in  one  of  the  existence  theorems  to  classes  U   of 
admissible  vectors  and   V  of  variational  vectors  as 
outlined  above,  we  prove; 

ThepreTTi  6:   Given  a  class  of  admissible  vectors   XJ  which 
contains  at  least  one  element,  and  given  a  fixed  element 
Vq  of  the  variational  class,  the  variational  function 
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P[U]  =  I  (U  -  Vq,  U  -  Vq) 

(8.2)  T  n  . 

=  I  ^  |U  -  VqT  dV 

is  minimized  on  the  class  of  admissible  U  by  U  =  Z. 

In  Theorem  6,  it  will  often  be  convenient  to 
replace  the  functional  by 

(8.3)  F[u3  =  I  (u,u)  -  (n,VQ)  . 

The  replacement  is  permissible  since  the  two  functlonals 

differ  only  by  the  constant  i  ("^0*^0^*   ^°^  ^^®  proof  of 

the  theorem  we  shall  require  the 

Lemma t   The  vector  U  -  Z  is  orthogonal  to   V  -  Z  for 

every  choice  of   U  and   V. 

The  vectors  U-Z   and  V-Z   are  easily  seen  to  satisfy 

homogeneous  conditions  and  the  orthogonality  requirements 

of  Table  II,  page  71,  are  satisfied  and  yield  the  lemma 

directly.   In  the  example  above,  for  Instance,  take 

X  =  U,   Y  =  V.   We  have 

curl  (U  -  Z)  =  0  ,      (U  -  Z)^  =  0 
div  (V  -  Z)  =  0  ,     [S;V-Z]  =  0 

so  that  orthogonality  is  direct  from  Table  Ila. 
We  then  prove  Theorem  6  as  follows? 

P[U]  -  P[Z]  =  I  (U-V,  U-V)  -  I  (Z-V,  Z-V) 


=  i  (U-Z,  U-Z)  -  (U-Z,  V-Z) 
=  i  (U-Z,  U-Z)  >  0 
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We  conclude  that  P[U]  >  P[Z]   and  the  equality  prevails 
only  when  U  =  Z. 

It  should  be  remarked  that  although  compatibility 
conditions  do  not  appear  explicitly  in  this  argument 
they  are  implicit  in  the  assiimptlon  that  at  least  one 
U  and  one  V  exist  satisfying  their  respective 
conditionso  As  a  matter  of  fact,  the  compatibility 
conditions  in  Theorems  1-3  are  characterized  abstractly 
by  the  simple  requirement  that  there  exist  at  least  one 
member  of  each,  the  curl  class  and  the  dlv  class* 
Conjugate  variational  problems o   The  Lagrange  multiplier  rule, 

Let  X  and  Y  now  denote  vectors  corresponding  to 
the  categories  (a)  and  (b)  of  Table  V  without  regard  to 
which  is  admissible  and  which  variational*   The  two 
functional s  for  the  conjugate  problems  are  given  in  the 
form  (80 3)  as 

P3^[X]  =  I  (X,X)  -  (X,Yq) 


(Soil) 


P2[Y]  =  i  (YsY)  -  (Y,Xq) 


for  X  and  Y  respectively,   Xq   and  Yq  being  fixed 
elements  of  the  respective  classeso   These  functionals 
have  a  common  minimizing  vector,   Zo   Adding  the  tw® 
minima  and  using  the  lemma  of  Theorem  6,  we  obtain 


(8,5)        P^^EZ]  +  P2[Z]  =  -  (Xq,Yq) 
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Since   (Xq,Yq)  is  fixed,  this  formula  tells  us  that  an 
upper  estimate  of  the  rninimum  value  in  one  problem  yields 
a  lower  estimate  for  the  conjugate  problem.   Pvirthermore, 
if  Xq  and  Yq  can  be  chosen  ©rthogonal  (say,  by  choosing 
one  to  be  identically  zero  if  it  satisfies  homogeneous 

conditions),  we  have  the  classical  theorem  in  which  one 

21 
miniimom  is  the  negative  of  its  conjugate   • 

As  a  simple  example,  let  us  take  a  case  from  the 

corollary  to  Theorem  1,  the  problem  with  closed  curve 

periods  J 

admissible t   div  U  =  0  ,  ^^  ~  ° 

(8.6)  ^   ^ 

variational!   curl  V  =  0,   [  nj^lVJ  =  v^     • 

Prom  (6»5)»  we  have  for  any  U  and  V  in  these  classes, 

(8.7)  (u,v)  =^  [rk;v]  [Zki^l 

so   that   the   functional    (8.3)   for  the   problem   takes   the 
form 

(8.8)        P[U]    =Jl^dY-^v^J     U  .  dS    . 
D  ^-^  Zk 

The   fixed  vector     Vq     no   longer  appears   in  the   fixnctional. 
It   is    represented  by  the   parameters      v^     which  deterraiae 
the   variational   class   explicitly  and   it    is  no   longer 
necessary  in  this    case  to   choose   a  specific  vector     Vq 
2l7      Courant-Hilbert    (loc»  clt.j.  ',  ~ 
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which  takes  the  values  of  the  periods  to  define  F[TJ],   If 
the  open  surface  representation  Yl     ^^   used  to  determine 
the  ///  group,  the  problem  then  takes  the  form 


admissible;   div  U  =  0  ,  U  =0,   [H^-jU]  =  J, 


2 
"^  dV 


(8.9)  {  variationals   curl  V  =  0 

Flu]  =  j  i  n' 

D 

Here  the  variational  function  is  simplified  by  taking 
as  the  chosen  element  of  the  variational  class,   V^-  =  0. 
Again  the  variational  function  does  not  depend  on  V^« 
Problem  (8.9)  is  similar  to  (8.6)  in  that  all  conditions 
are  homogeneous  except  for  a  set  of  //C -periods,   but  in 
one  case  these  periods  have  been  specified  as  closed  and 
in  the  other  as  open.   This,  of  course,  necessitates 
transferring  the  period  conditions  from  the  curl  class 
to  the  div  class  in  order  to  keep  the  grouping  of  Table  V 
inviolate.   Now,  applying  the  classical  procedure  of  the 
method  of  Lagrange  multipliers,  we  alter  problem  (8,9) 
and  look  upon  div  U  =  0,  U^  =  0,   as  admissibility 
conditions  while   [5~i,jU]  =  ii,   is  a  side  condition 
(this  is  neither  an  admissible  or  a  variational  condi- 
tion).  Now,  using  the  Lagrange  multiplier  rule  we  obtain 

(8.10)  P[U]  =  f  I  U^  dV  -  ^  -K^      J      U  .  dS 


k=l 


IL 


k 
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as  the  variational  function.  The  two  problems  (i.e., 
with  open  and  closed  periods  specified)  have  been  put 
in  the  same  form,  but  the  functlonals  (8.8)  and  (8.10) 
have  slightly  different  interpretations.   In  (8.8)  the 
constants  v^     which  appear  in  P[U]   are  given  data. 
In  the  Lagrange  multiplier, problem,  the  constants  X, 
which  appear  in  the  functional  (8.10)  are  not  prescribed 
but  must  be  determined  afterwards  by  the  side  conditions 
so  as  to  yield  the  given  values  of  X  . 

The  example  is  directly  applicable  to  a  magneto - 
static  field  U  =  B.   If  the  boimdaries  are  perfect 
conductors,  with  B   =  0   as  the  appropriate  boundary 
condition,  the  period   [  P',  iU]   represents  an  m.m.f.  or 
the  current   j,   in  the  linked  circuit   PT.  while 
[J",  }U]   is  simply  the  flux  linking  the  circuit  loop 
P,  .   If  the  fluxes  are  prescribed,   [J^,  ]  =  J^,   then 
B  minimizes  the  magnetic  energy  J  Pu  ^   ^^*   "^"^  ^^^ 
currents  are  given,   [  P,  ]  =  P-Ji,*   then  B  minimizes 

the  functional  J   ^  B^dV  -  i±  YZ.   Jir  Iir^B^*   ^^  *^®  °^® 

D  ^  k=l  ^  ^ 

case,  admissible  vectors  B  must  be  chosen  consistent 

with  prescribed  values  X  ^   -^^  ^^®  other,  the  given 

values   1   do  not  affect  the  admissible  class  but 

appear  as  natural  boundary  conditions  on  the  minimizing 

vector.   The  third  possibility  (the  Lagrange  multiplier 

case)  is  to  consider  the  fluxes  as  given,  but  to  ignore 
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this  restriction  in  choosing  admissible  vectors,  use 
the  second  variational  function,  and  adjust  the  7^., 
to  obtain  the  given  X  • 

We  have  seen  in  this  example  how  the  replacement 
of  an  admissibility  condition  by  the  dual  condition 
(according  to  Table  V)  on  the  variational  side  leads 
to  a  change  in  the  functional  of  the  problem  describ- 
able  in  terms  of  the  Lagrange  multiplier  rule.  This 
is  frequently  the  case,  but  sometimes  requires  a 
generalization  of  the  Lagrange  multiplier  rule.   For 
this  reason,  the  formulation  of  Theorem  6  may  be 
considered  to  generalize  as  well  as  justify  the  Lagrange 
multiplier  rule.   (See  the  appendix). 

Variational  formulation  of  the  problem   curl  curl  A  —  -^^ 
dlv  A  =  cT  , 

Let  us  consider  one  more  special  problem. 

I^dmissiblej   curl  U  =  j^. 

(8.11)/ 

/variational;   div  V  =  0,   V  given,   [2_}V]  given. 

Let  us  prescribe  in  addition  that  V   satisfies 

n 

[S,  JV]  =0   so  that  a  regular  vector  potential  with 
V  =  curl  A  exists.   From  (6.18)  we  have 

(U,V)  =  (A,;;^)  +  J    A^x  U^  •  dS 

a* 

whence 
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(8.12)  P[TJ]    =  Jl  u2  dV  -  /    A^kU^  .dS 


where   the    term      (A,   ]^)     being  constant,    is  dropped. 
The   inforraatlon     V       given  on     S       n©w  appears   in  the 
functional   disguised  in  the    form  of   the  multiplier     A. J 
where  we  note   that     A^     is  not  uniquely  defined  by  the 
given     V   ,     but   that   any  representation  suffices. 

By  taking  the   conjugate  problem,   but   keeping   the 
period  conditon  ©n  the   variational   side,  we   obtain  a 
more    Interesting  situation, 

admissible t      div  U  =  0,     U       given 

(8.13)  1  1  r^ 

variational!   curl  V  =  j^,   [  KJV]   given. 

Under  the  prescription  on  U  ,   IS,  JU]  =  0,  we  have 

U  =  curl  A,   The  functional  appropriate  to  this  problem, 

when  expressed  in  terms  of  A,   takes  the  form 

(8.14)  F[A]  =  J||curl  a|^  dV  -  f  A^   •  IE  dV  -  y  A^  ><  U^  •  dS  , 

There  exists  a  unique  minimizing  vector  U,  but  the 
vector  potential  A  is  not  unique  since  an  arbitrary 
gradient  may  be  added.   To  assure  the  uniqueness  of  A, 
we  may  proceed  in  several  ways  by  referring  to  the 
loniqueness  theorems.   Supplementary  conditions  may  be 
imposed,  subject  to  compatibility  requirements,  such  as 
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(8,l5)        div  A  =  <r  ,  A^,,[SJA]  or  [C|A]  given 

The  problem  will  take  the  desired  form  when  the 
^-periods  are  given  in  the  form   [ZIJ^^   rather  than 
[PjV],   thereby  putting  the  period  condition  on  the 
admissible  side.  The  boundary  values  A.   must  then 

A 

be  chosen   compatible  with     U        and  with  periods      [  P'jA] 
compatible   with      i^iU]*      Instead  of   requiring   this 
compatibility  condition,   we  may   delete    all  reference 
to     U     from  the   problem  and  leave   the    information  in 
A,      whence     U       and      [^{U]      may  be  computed.      This 
yields   a  problem   in  which     A.       can  be    arbitrarily 
specifiedl 

fadmissiblef      div  A  =  gT,   A^.,    [SjA]    or   [CjA]  given 
(8.l6)f  ^  _ 

[variational J      curl   curl  A  =    — 

with   the    variational   fixnction  given  by    (8.11}.).      The 
surface    integral   Is  now   a  constant   of   the   problem  and 
may  be    deleted,    yielding 

(8.17)  P[A]    =  /I    Icurl  A|^  dV  -   /a  '^^  dV  . 

D  D 

In  the  problem  (8.16)  -  (8»17)  we  have  found  a  vari- 
ational formulation  of  Theorem  li.« 

We  remark  that  the  variational  problems  (8.11)  - 
(8,12)  and  (8,16)  -  (8,17)  are  conjugate  (in  a  somewhat 
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generalized  sense,  sinoe  one  involves  U  and  the  other, 
the  vector  potential  ©f  V)   since 


Pn tU]   +  Pp^^^  ~  °®^®** 


whe 


re  P^  is  defined  by  (8 ,12)  and  Pg  by  (8.16), 
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9e  A  General  Existence  Theorem  and 

22 
Related  Considerations   « 

The  existence  theorems  1-3  may  be  extended  to 

mixed  boimdary  value  problems  In  which  either  one  of 

Z„  or  Z.   Is  prescribed  on  each  of  the  disjoint 
n       u 

connected  components  S,   of  S   but  not  necessarily 

over  all  S  ,  We  shall  denote  by  S   that  section  of 

S   on  which  normal  conditions  Z   are  specified,  and 

n 

by  S.   the  remainder  of  S  ,   and  there  we  specify 

tangential  boundary  conditions  Z^.   In  particular.  If 

S   Is  connected,  either  S   or  S.   is  vacuous,  and 

n       u 

we  have  one  of  the  previously  considered  cases, 

it  *k 

In  the  event  that  neither   S   or  S.   is  vacuous 

n       \0 

it  is  no  longer  true  generally  that  the  open  periods 
are  determined  on  a  homology  class  by  the  value  on  a 
single  member.   For  example,  in  Figure  6,  page  33,  let 
St  =  S  ,  So  =  So,,   If  >  T   is  a  member  of  the  class 
of  '5~-i   such  that  the  two  surfaces  can  be  completed 
on  S   to  a  closed  bounding  surface  in  D  only  by 
the  addition  of  part  of  Sp  (i.e.  of  S^^),   then  it  is 
not  possible  to  specify  the  integral   [^^;Z]   in  terras 

of  the  integral   [y~-.  jZ]   since   Z   is  not  specified 

^  t 

on  Sp»   Similarly,  for  open  arcs   C   and  C   connecting 

22*   This  section  may  be  omitted  or  postponed  without 
disturbing  the  continuity  of  exposition. 
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the  two  botmdary  components  but  initiating  at  different 
points  on  S^  the  integral  [0  JZ]  cannot  be  expressed 
in  terms  of  the  integral   [CjZ]   since  X^  is  not  given 
on  S  •  Nonetheless,  by  a  simple  device  it  is  possible 
to  recover  the  property  that  the  periods  on  a  homology 
class  can  be  obtained  from  the  period  on  any  representa- 
tive ©f  the  class  and  the  known  bovuadary  values.  The 
method  is  to  utilize  a  mixed  representation  for  the 
abstract  group  consisting  of  both  open  and  closed 
elements  and  to  extend  the  concept  of  homology  modulo 
the  boundary. 

Let  us  first  consider  the  /ri    group  in  its  dual 
representations,  \  f"' j  the  group  of  closed  curves  in 
D,   and  )2Ij  ^^^    group  ©f  ©pen  sxirfaces   [mod  S  ]. 
Let   ["*■]''•••>  n'   denote  a  basis  for  \V^j   which  we 
may  assume  to  consist  of  simple  closed  curves  on  S 
and  let  T"t  >•«  vT"   be  a  reciprocal  basis.   The  most 
easily  comprehended  method  of  forming  a  mixed  basis  is 
to  take  a  subset  of  the  X^*f      say  \^^j»»*    ^]^'      ^"-^ 
t©  complete  the  basis  with  ^Zir+i^*  ** '2^*  ^   ^®  given 
boTindary  value  problem   1"^,  ;,•••.  PV.  are  chosen  as  a 
basis  of  that  subgroup  of  |  \^l     which  is  spanned  by 
the  closed  ctirves  of  S^«   The  abstract  subgroup  for 
which   r^-,,,,.rT   is  a  closed  curve  basis,  is  denoted 
by  /vf.  •   The  orthogonal  complement  t©  ^^»      which  we 
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denote  by    7?^.  ,      is   generated  by   a  basis      l"^k+l^  ••• -"  ^t) 

'Sf' 

which  lies  completely  on  S   since  we  have  assiamed  each 

basis  element  to  be  connected  ^«  The  other  representation 

of  7yi±.     in  terms  of  the  basis  of  ©pen  surfaces 

y,  J,  ,...  T"   l^as  also  a  simple  interpretation  in  terms 
•^ — K+l      ' — ^p 

of  homology.   It  is  intuitively  seen  that  the 

ZIb-+i  ?  •  •  •  ?XIt)  "**y  ^®  chosen  so  as  not  to.  intersect 
S. •  The  group  ^.   may  therefore  be  redefined  as  the 
homology  group  of  ©pen  stirfaces  in  D  modulo  3  J 
that  is,  the  group  of  open  surfaces  with  botindarles  •n 
S  •   Open  surface  periods  are  now  defined  on  }/l. 
exactly  as  before.   Given  an  element  in  a  homology  class 
of  iTtft      the  open  surface  periods  for  all  other  members 
can  now  be  obtained  from  the  known  boundary  values  on 
S  .   For  example,  in  Figure  ?♦  page  36,  let  S   consist 
of  S^  and  S,  •  The  curve   PL  may  be  taken  as  the 
basis  of   ^^>   the  s\arf  ace  5"..   bounded  by   PL   as 
the  basis  for  ^K^* 

Similarly,  for  the    G  group  we  let  C^     consist 
of  the  subgroup  spanned  by  the  closed  surfaces  constituting 

24*  Note  that  the  orthogonal  subgroup  to  T^tZ  may  be 
smaller  than  the  subgroup  which  is  spanned  by  all 
closed  curves  of   f Ij   on  S*«   For  example,  in 

Figure  6  with  S*  =  S*,  S^   =  S'^,  M^     has  the 

basis  ^pt    r^p  ^hil®  the  orthogonal  complement 
is  the  null  apace* 
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it  -5— 

S  •  The  ortiiogonal  complement  ^   may  be  assigned  a 

...  )  ^ 

basis  consisting  of  all  the  components  of  S.   but  on©« 

The  group   £,   may  also  be  thought  of  in  its  other 
representation  as  the  homology  group  ®f  open  curves 
[mod  S.  J   with  a  basis  consisting  of  open  arcs,  each 
connecting  the  excluded  component  to  one  ©f  the  other 
components  of  S.  . 

We  are  now  able  to  state  the  proper  generalization 
of  Theorems  l-3e 
Theorem  7;   Let  D  be  a  domain  with  boundary 

S*^  =  S*  +  si!  with  S*   and   st  disconnected.   There 
n    t         n        t 

exists  a  unique  solution  of  (5ol)  in  D  taking  boundary 
values  X   on  S   and  X^   on  S^   compatible  with 
cr   and  ^,   but  otherwise  arbitrary,  and  with  arbi- 
trary values  for  a  set  of  periods  on  7?t^     and  ^  ^, 
Note  that  the  given  boundary  values  explicitly 
specify   [  P;Z]   on  '^^   and   [SjZ]   on  C^     so  'that 
complete  sets  of  periods  on  both   /H  and   G  are 
given  by  the  data. 

If  neither  S   or  S.   is  vacuous  then  there  are 
n       "C 

no  compatibility  restrictions  on  Z^.   The  requirement 

—        ■» 
on  Z,   is  that  it  be  compatible  with   -.   °^  ^f 

Specifically  if   P  bounds  21*  °^  ^   ^®  ^^^® 

r  P»Zl  =  J^    ~^     dS.   Here,  as  before  in  the  inhomogeneous 

case,  it  is  necessary  to  utilize  generalized  periods 
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which  take    Int^j   accoimt   the   prescribed  boundary  values 
and  the   prescribed  fields     fT      and     "^t 

Next,   we  derive   the  orthogonality  relation  for  the 
system    (6.2)   with  the  mixed  boundary  conditions 

X,  =0     on     S,  =gS* 

(9.1) 

Y^  =  0     on     S*  =  J±Z  Sj 
^  "       k=v+l     ^ 

Observing  that  ^     is  constant  on  each  component  of 
S,   we  obtain 

(9.2)  (X,Y)  =fZ^i,   i\i^^    +^  tr^ki^^  EkJ^^ 

k=0  k=l 

By  writing   (9.2)   in  terms   of  the  homology  groups     <£, 
and     7?^       we   obtain   an  Invariant  form  generalizing 
Hodge's    theorem, 

(9.3)  (X,Y)   =11   [G^iXl    [S^JY]    +^   [r^JX]    CZj^jY]    , 

^n  ^t 

in  which   the   summations   are   taken  over  basis   elements 
of   the    Indicated  groups.      In  other  words,    the    inner 
product  of  Irrotatlonal     X     and     solenoldal     Y     over     D 
subject    to    (9.1)    is   the   bilinear  product  of   the    two 
types   of     ^.      periods   plus    the  bilinear  product   of 
the   two    types   of      ^         periods. 

Vectors     X,Y     satisfying   (6.2)   and    (9.1)  will  be 
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orthogonal  if 


(9.1^) 


[C^iX]   =  0  or   [S*;Y]  =0  on  ^^ 
[  Pj^jX]  =  0  or  CZkJYJ  =0  ®^    Kt 


With  the  aid  of  these  orthogonality  conditions  it  is 
possible  to  further  refine  the  decomposition  of  harmonic 
vectors  given  in  Table  III,#  The  most  general  such 
decomposition  is  the  following; 

Table  VI 


1 
H 

ti 

H 

t                        ^ 

H.    =  0  on  S^ 

[  PjH    ]    =  0  on  }^^ 
[C;h']    =  0  on  "g:^ 

H     =  0  on  S„ 
n                    n 

[XI;  H  ]  =  0  on  7rl^ 
[SJH   ]    =  0  on   ^^ 

A  decomposition  of  this  type  exists  for  each  separation 
of  S*  into  st  and  S^  and  for  each  further 

decomposition  7^^  ""  '^t  "^  ^  t  ^^  6n  "  ^n  "^  ^n' 
That   h'   and  H   are  orthogonal  follows  immediately 
from  (9,i|.).   Given  an  arbitrary  harmonic  H,  Theorem  7 
yields  a  unique   h'   satisfying  the  criteria  of  Table  VI 
and  the  conditions   H^  =  H^  on  S^,   [ZIJH  1  =  IZ-*'^^ 
on  yjfl      and   [Sjh' ]  =  [SjH]   on  "JJ^.   The  vector  H 
is  given  by  subtractions* 
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The  decomposition  given  by  Table  VI  can  be  described 
as  follows.   Given  certain  data  which  define  a  unique 
harmonic  vector  according  to  Theorem  7,  we  split  the 
boundary  data  and  each  type  of  period  data  into  two 
sets  such  that,  when  supplemented  by  complementary 
homogeneous  data,  each  set  determines  a  unique  harmonic 
vector,  and  the  two  are  orthogonal. 

To  prove  Theorem  6  in  general  requires  only 
slightly  more  effort  than  in  the  special  case.  Consider, 
for  example,  the  variational  problem  in  which  the  data 
are  divided  as  follows  without  specific  reference  as  t© 
which  class  is  admissible  and  which  is  variational. 

curl  X  =  _2L,  X^  given  on  S^,  [  HjX]  given  ©n 

TKI*    tC|X]  given  on  "Jf^ 

(9.5)  / 

\  dlv 


V 


Y  =  cT  ,   Y     given  on  S    ,    i^iY]    given  on 
>^"'    f^*^^    Slven  on  "^". 


To   exhibit    the   orthogonality  of     X— Z      and     Y-Z     we 
obtain  from  Tables   IV  and  VI 

(9.6)  Z  =  curl  A+H     +H     +gradjZ) 


and 
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X  =  c\irl  A  +  H  +  H  +  grad  )8 


I   -» 


Y  =  curl  A  +  H  +  H  +  grad  0 

where  the  unbarred  symbols  represent  quantities  which 
are  uniquely  determined  by  the  data  in  (9»5)  defining 
the  classes  X  and  Y,   and  the  barred  symbols  depend 
upon  the  particular  choice  6f     X  and  Y»  Per  example, 
expanding  Y  according  to  Table  IV,  the  first  column 
alone  takes  the  condition  dlv  Y  =  ar-  and  this  with 
the  other  conditions  in  the  column  determine  grad  f6 
uniquely.  Prom  the  decomposition  (9.7)  It  is  clear 
that  Z  is  the  only  element  common  to  both,  classes 
X  and  Y.   The  orthogonality  (X-Z,  Y-Z)  =  0  is 
immediate  from  (9o6)  and  (9.7).  Theorem  6  Is  thereby- 
proved  for  the  general  Theorem  7.   The  problem  with  mixed 
boundary  conditions  can  now  be  put  as  a  variational 
problem  of  the  type  of  Section  8» 
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10,   The  Electromagnetic  Equations. 

We  shall  now  apply  the  existence  theorems  of 
Section  5  to  the  formulation  of  well-posed  problems 
for  the  pre-Maxwell  electromagnetic  equations.   First 
by  employing  c)B/()t  instead  of  B   as  the  fundamental 
variable,  the  application  of  the  existence  theorems 
is  made  immediate.   Explicit  time  dependence  is  next 
introduced  by  means  of  Paraday' s  law.   Most  significant 
for  fluid  magnetics  is  the  time  dependent  problem  in 
which  only  enough  data  are  prescribed  to  determine   B« 
The  results  are  then  applied  to  a  few  simple  examples, 
namely  Faraday's  law  for  a  physical  wire,  the  homopolar 
generator  and  the  pinch  effect. 
Well-posed  problems  for  E  and  B. 

We  shall  not  begin  directly  with  the  pre-Maxwell 
system  (1.1)  but  instead  take  the  simpler  system  obtained 
by  removing  time  from  the  equations,  taking  B  =  jB/^t 
as  the  fundamental  variable  instead  of  B.   Differ- 
entiating with  respect  to  time  in   div  B  =  0   and 
curl  B  =  ^lJ,   we  obtain  the  system 


curl  E  =  -B       curl  B  =  uJ 
(10.1)      <> 

div  E  =  q/|c       div  B  =  0 


where   J  and   q   are  assumed  to  be  given.   The  time 
dependence  has  been  eliminated  in  the  sense  that  we 
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can  solve  for  B  and  E  at  any  fixed  Instant,  given 

appropriate  boundary  conditions  and  period  data,  provided 

that  at  the  given  time   q  and  J  are  knovm  as  functions 

of  space.   For  example,  boundai^  conditions  on  B  might 

be  obtained  by  differentiating  prescribed  boundary  data 

for  Bo   The  situation  here  is  similar  to  that  In 

Irrotational  incompressible  fluid  flow  where  the  flow 

is  determined  at  each  Instant  by  the  instantaneous 

boundary  values  and  varies  in  time  only  if  the  boundary 

conditions  change.   An  interesting  special  case  in  which 

this  formulation  is  identical  with  the  system  (1.1)  is 

that  for  which  time  may  be  separated  out  either  sinus- 

oldally  or  exponential ly. 

« 
With  the  observation  that   curl  curl  E  =  -|xJ  we 

have  from  Theorem  k-i 

Theorem  8f   There  exists  a  unique  solution  of  (10.1 ) 

for  E   and   B  possessing  arbitrary  boundary  values 

E.   and   C-perlods  in  either  representation,  closed  surface 

[SjE]   (compatible  with   q)  .^or.open  curve   [C}Ej^ 

Since   -E   is  the  vector  potential  of  B,   the 

problem  (10,1)  is  the  same  as  the  variational  problem 

of  (8,11)  et,  seq,  for  B«   In  particular.  Theorem  8 

corresponds  to  the  variational  problem  (8,16)  -  (8.17  )• 

To  recapitulate,  the  solution  for  E   is  given  by 

minimizing 
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(10.2a)  P[E]    =  J  i    Icurl  EJ^  dV  +  J   E«  tiJ  dV 

D  D 

subject   to    the  admissibility  condition 

(10.2b)   dlv  E  =  q/Ki,   E.    given  on  S   ,    [S}E]    or    [GjE]    given, 


The  value   of     B      is    then  given   as     B   =   -curl  E,      or   it 
may  be   found  directly  as    the    solution  of   a  problem  of 
the  form   (8.11)   -    (8.12)   by  minimizing 


(10.3a)  P[B]    =  J  i  B^  dV  +y    E.  '^  B  •  dS 

D  ^  S, 

subject   to    the   admissibility  condition 

(10.3b)  curl  B  =  p.J  • 

Only  the  given  boundai^r  values     Ex^      appear   in  the 
variational   function     P[b3      which   is    therefore  well- 
defined. 

A   second   existence    theorem   is   obtained  by  applying 
Theorem  3   to      B      and     E      in   turn. 

Theorem  9?      A  unique    solution  of    (10.1)   exists  with 
given  boiindary   values      B        and     /v^-periods      (CXI*^^      °^ 
[  n'jB]  )      as  well   as   given  boundary  values     E        and 
^-periods      ([XI;E]      or      [PjE]). 

In   physical   applications    some   of   the    period 
conditions    are    seldom   encoimtered.      For   the   magnetic 
field  the     /^-periods  in  both, versions,    the  m»mef»'s    s 
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[  PjB]      and  fluxes      [^Z*^^      ^^®   generally  relevant  to 
applications,   while   the      C -periods   either  vanish  in 
the  form     [S;B]      or  are   seldom  relevant   in  the   form 
[CjB]      as   difference   in  magnetostatic   potential.      For 
the  electric   field  the  e«m»f»'s      [P'JE]      and  charges 
[S;E]      are  both  generally  relevant,   while   the  open  periods 
[C;E]      as  differences  in  electrostatic   potential   are 
sometimes  relevant   and  open  surface  periods      [^"fE]      as 
electric   fluxes  apparently  never  are* 
Well-posed  problems   for     E     and     B« 

Returning   to   the   electromagnetic   equations    (1»1) 
with     B     as  variable  rather  than     B     we  consider  the 
crucial   equation 

4|  +  curl  E  =  0 

as  being  the  only  one  in  which  time  variation  enters 
explicitly  and  also  the  only  one  in  which  E   and  B 
are  coupled.   Let  ^  t»e  a  surface  moving  at  velocity 
u  and  let  [~'     be  its  boiinding  curve.   Prom  the  Identity 

(lO.i^.)  ^  /  Z  •  dS  =  J     (^  +  u  div  Z)  .  dS  -  y"  uxZ  •  dx 

for  the  rate  of  change  of  the   flux  of   an  arbitrary 
vector  field     Z     through     S,      setting     Z  =  B,     we   obtain 
Faraday's   law, 
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(10.5)^^     -l^    y  B  •  dS  =  -  /  (E  +  u  »B)  •  dx  . 
dt  ^  p 

We  note  that  previously  the  analysis  held  as  well  for 

moving  boundary  as  for  fixed  boundary J  any  motion  of 

the  boundary  was  irrelevant;  the  analysis  applied  wholly 

to  the  configuration  which  obtained  at  a  given  instant. 

With  the  utilization  of  (10«5)  on  the  boimdary  of  the 

domain,  however,  compatible  boundary  values  of  B  and 

B  involve  the  motion  of  the  boundary.   Nothing  new  is 

obtained  from  (10.5)  when  applied  t©  the  interior  of 

the  domain;  it  merely  expresses  a  property  of  any  solution 

of  the  equations. 

P©r  a  closed  curve   f^  on  S   there  are  two  cases 

to  consider!   either  the  surface  ^  bounded  by   P*  may 

be  taken  as  part  of  S    (that  is  \^     is  homologous  to 

zero  on  S  )   or  ^Z  belongs  to  one  of  the  homology 

classes  of  ^T"  |  .   If   P  '^  0  on  S^,   (10.5)  relates 

the  boixndary  values   (E  +  uxB).   to  the  values  of  B  • 

"C  n 

. 

In  fact,   B   is  completely  determined  by   (E  +  uxB)., 

. 
but  if  B    is  known  as  a  function  of  time,   (E  +  u>^B)+. 

is  deternined  only  to  within  an  arbitrary  additive 

surface  gradient.   In  the  second  case  with  ^  in   <^/> 

25.   An  extension  of  (10.5)  to  open  arcs  which  is 

sometimes  useful  is  the  following  easily  derived 
formula 

^  f  {A*  Ax)  =  -  f   (E  +u'^B)  -dx  +J'd(u»A). 
'^^  C  C  C 
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we  conclude  that  boundary  values  of   (E  +  uxB).   are 
related  to  the  values  of  the  periods   [^|Bl,   specifically 

^  [HkJB]    =  -    ir^iE  +U.B]    , 

/\ 

where      p,       is    the  boundary  of     y~.       in  the  notation 

introduced  in  Section  3«      In  other  words,    if     B^      and 
[^JB]      are  known  at     t  =  0     and      (E  +  u<B)^     is  known 
for     t  >  0,      then  we   can  compute      B        and      [^IjB]      f©r 
t  >   0»      The   converse   is  not   true.      Knowledge   of     B       for 
all   time   only  yields      (E  +  u^B)^      to  within   an    arbitrary 
surface   gradient    and      [JIjB]      fixes   only  h:.lf  of  the 
periods   of  this    surface  gradient,   namely,      iV^]     but 
not      [  P] . 

It   is    sometimes  helpful  to    adopt   the    interpretation 
of      (E  +  u^B),       as    the    surface   e.m.f.   measured   in  a 
coordinate    system  which  is    instantaneously  moving  with 
the   local   velocity  of  the    surface   at    the  point   observed. 
We   recall   that   this   is   a  Galilean  not   a  Lorentz   invariant 
interpretation   (cf.    footnote ,2, , p.    10)# 

The    following  immediate  modifications   of  Theorems   8 
and   9   apply  directly  to    the   pre-Maxwell  equations    (l,l)« 
The  boundary,      S    ,      of   the   domain   is   assumed   to   have   an 
arbitrary  but   given  velocity     u     which  can  vary  over    the 
surface   and  in  time. 
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Theorem  10?      A  unique   solution  of   the   pre-Maxwell 

system   (1.1)   for     E     and     B     for     t  >  0     exists, 

taking    given  boundary  values     E^      or      (E  +  u^B). 

and  periods      [SjE]      (or      [C;E]      in  case   the  botmdary 

condition   is     E^^ ) ,      all   assigned   arbitrarily  for      t  >  0, 

together  with  boundary  values     B       and     //(^-pepiods, 

[PjB]      or     [2I;B],      assigned  only  at      t   =  0.      The  only 

compatibility  condition  is    J^  B     dS   =  0   • 

S   " 

If  E.   is  given,  we  refer  to  Theorem  8  to  solve 

for  B  and  E,   and  then  compute  B  from  its  value 
at   t  =  0  which  is  determined  by  the  given  initial  data 
(Theorem  3).   If   (E  +  ux-B),   is  given,  we  first  compute 
from  the  known  initial  values  the  boundary  values   B 
and  periods   [X!*^^   ^^  functions  of  time,  as  discussed 
in  connection  with  (8.5)«   Then,   B  is  computed  at  each 
instant  using  these  values  and  Theorem  3«   Finally,   E. 
is  computed  as   (E  +  u  ''■B).  -  (uxB).   and  then  E  is 
obtained  by  Theorem  \\.,      Compatibility  is  trivially 
verified  at  each  stage. 

Theorem  11?   The  identical  statement  is  made  for  (1.1) 
as  in  Theorem  9,  with  B   replaced  by  B,   and  all 
boundary  and  period  data  specified  as  functions  of  time, 
t  >  0. 

Note  that  the  data  specified  here  exceeds  that 
specified  in  Theorem  9  by  the  inclusion  of  B    and  its 
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periods   at      t   =  0. 

Theorem  12>      A  \mique   solution  of    (1.1)   exists   for   the 
vector     B     If  we   specify     B        for     t  >  0,      [P|B]      or 
tZZi^^      at      t  =  0,      and   the  e.in.f.'s      [P^jE+uxB], 
(k  =  l,.,.,p)      for   all   time. 

This    Is    a  special  case   of  Theorem  10   In  which  only 
enough  data  has  been   specified   to   determine     B}      E     will 
not  be  unique.      In  particular,    we   conclude    that   Insofar 
as     B      Is   concerned,    the    surface   distribution  of   the 
"applied  voltage"      (E  +  uxB).       on  the  boundary   is 
immaterial;    only  total   e.m«f,'s    (and  only  half  of   these) 
are    relevant. 

An  important   general   remark   is    that,    despite   the 
appearance   of   a   time   derivative      ^B/(3t      in  the    differ- 
ential  equations,    one   cannot   give    initial  values   of     B, 
or  of    any  other  quantity   as   functions   of   space   over   the 
domain     D,      but   one   can  only  give   initial  values   of 
certain   constants    (periods)    and,    on  occasion,   boundary 
values. 

This    theory  will  be   Illustrated  by   a  few   simple 
examples.      Suppose     E.       is    a   given   surface   gradient   on 
the    surface   of   a  torus.      If     B      =  0      initially,   >ie 
conclude   that      B      =0      always.      The   complete   electro- 
magnetic field   inside    (outside)    the    torus  becomes 
uniquely  determined  when  we   specify  the   initial  value 
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of  the   total  magnetic   flirx  inside    (outside).      Of  the 

A 

two  given  periods   [  Pl^]   and  [PjE],   one  determines 

the  rate  of  growth  of  the  magnetic  field  inside  the 

torus,  and  the  other,  outside.  We  may  replace  the 

given  E.   by  any  other  surface  gradient  with  the  same 
t 

periods  without  changing  B|   the  electric  field  E 
in  D  will  change,  of  course.   In  other  words,  B  is 
insensitive  to  the  precise  way  in  which  E  is  applied. 

If  the  torus  is  in  motion  and   (E  +  uxB).   is  a 
specified  surface  gradient,  the  situation  is  exactly 
the  same.   In  particular,  if  the  e.m.f.'s  vanish. 


A 


[  P,  JE  +  uxB]  =  [  Pi,}E  +  u  "^B]  =  0,   then  the  two  fluxes 
(inside  and  outside)  are  constant  In  time.   We  can 
compute  B   at  each  instant  knowing  only  the  constant 
value  of  the  flux  and  the  instantaneous  shape  of  the 
torus;  the  velocity  of  the  torus  is  irrelevant. 

This  conclusion  is,  of  course,  not  restricted  to 
a  torus.   In  any  topology.  If  the  particularly  simple 
boxmdary  condition   (E  +  uxB).  =0   is  satisfied,  the 
magnetic  field  is  determined  at  each  instant  by  the 
instantaneous  geometrical  configuration  alone.   The 
electric  field  will  depend  also  on  the  instantaneous 
velocity  of  the  boundary. 
Perfect  conductors*   Faraday's  law  for  a  wire. 

The  property  of  vanishing  e.m.f. 's  is  attributed 
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to    a  perfect  conductor.      In  fact,   we   define   a  perfect 
conductor  precisely  by  the    constraint 

(10.6)  E  +  u  xB  =  0      . 

For  the  present  purpose,  we  assume  that  the  flow  velocity 
u(x,t)  of  the  perfectly  conducting  fluid  is  given 
a  priori}  in  fluid  magnetics  problems  it  will  be  determined 
together  with  the  electromagnetic  field  by  (1,1),  (10.6), 
and  appropriate  hydro dynamical  equations  of  the  fluid. 
Combining  (1.1)  and  (10.6),  we  have 

(10.7)  ll  =  curl  (u  xB)   , 

from  which  it  is  possible  to  conclude  that  B   is 
determined  uniquely  for   t  >  0  when  its  value  at   t  =  0 
is  given.   We  can  now  compute  E  from  (10.6)  and  so 
derive  the  values  of   q  as  well  as   J.   In  other  words, 
a  well-posed  problem  for  a  perfectly  conducting  fluid 
with  flow  velocity  known  consists  In  giving  only  the 
initial  values  of  B   and  requires  no  Information  about 
q  or   J.   In  particular,  the  values  of   q  and   J 
which  result  may  not  satisfy  conservation  of  charge 
(cf.  remarks  of  Section  l)« 

The  integrated  form  (10.5)  of  d^/^t   +  curl  E  =  0 
states  that   (E  +  u  kB).   is  continuous  across  a 
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26 
discontinuity  surface   •   In  particular,  it  is  zero  next 

to  a  perfect  conductor,  and  this  can  be  taken  as  an 

appropriate  boundary  condition  for  a  domain  (e.g., 

vacuum)  bounded  by  a  perfect  conductor.   Prom  the 

integral  form  of  dlv  B  =  0  we  conclude  that  B   is 

n 

continuous  across  an  interface,  and,  therefore,  that 

the  magnetic  flux  through  an  element  of  surface  is 

constant  in  time  at  a  given  point  on  the  perfectly 

conducting  boundary,  i.e.,  at  a  point  moving  with  the 

boundary  (cf,  (10.?)).   However,  this  boundary  condition 

is  not  independent  of  the  previous  one  and  could  be 

deduced  immediately  from   (E  +  uxB),  =  0   and  (10#5) 

as  the  compatible  boundary  condition  on  B  .   In 

particular,  the  compatible  boundary  conditions 

B  =  0,   (E  +  u  <B),  =  0   are  frequently  used  at  the 

surface  of  a  perfect  conductor. 

A  moving  perfect  conductor  provides  a  means  of 

actually  constructing  a  physical  situation  which  is 

appropriate  to  Theorem  12  in  all  generality.   If  we 

have  a  perfectly  conducting  fluid  (e.g.  a  dish  of  mercury) 

with  B   j^  0,   then,  insofar  as  the  motion  can  be  con- 
n     ' 

trolled,   E,  =  -(ux>B).   can  be  given  quite  arbitrai*y 

values, 

26.   We  are  accepting  the  integral  form,  as  basic  rather 
than  the  differential  form. 
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Let  us  now  consider   an  illustration  of  Faraday's 
law,    (10,5)   in  the   domain  exterior  to    a  conducting  wire. 
As   stated,    (10.5)  refers   to   the   flux  and  e.m.f.    of  a 
inathematical   closed  curve.      In  practice,    the   law  is 
frequently  used  for   a  physical   wire  which   is    almost   a 
closed  curve,    but   is    actually  cut    at  one   point  for   the 
insertion  of  a   voltage-measuring  device.      For   simplicity, 
ass\ime   the  wire   is   a  thin   (but  not   one -dimensional) 
perfect   conductor,    and    the    arc    joining  the   ends  of   the 
wire,      A     and     B,      (Figure   12)   is  not  moving.      The  wire 


Figure   12 »      Conducting  wire   loop, 

cannot  be  made   one-dlmenslonal   since    the   differential 
equations    (1,1)    do   not   allow  boiindary  values    to   be 
imposed  on  a  one-dlmenslonal  manifold.      This  has   the 
desirable   consequence  .that   the   perturbation  of   a  given 
electromagnetic    field  due   to   the   insertion  of   a  thin 
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wire  Is,  roughly  speaking,  concentrated  within  a  few 
diameters  of  the  wire.   Specifically,   J  =  t/  B  •  dS 
is  almost  the  same  for  the  original  and  for  the  perturbed 
field.   Since   (E  +  u  xB)  =  0  within  the  perfectly  con- 


ducting wire,  we  have 

B 


/  E  •  ds  =  ~ 


dj 

dt 


where  the  integral  is  taken  along  the  short  arc  AB 
which  closes  the  circuit.   In  other  words,  although  the 
electromagnetic  field  Is  perturbed  only  near  the  wire, 
it  is  very  radically  changed  there,  and  although 
J^(E+U'<B)«dx  Is  hardly  altered,  the  distribution 
along   P  of   (E  +  u<B).   is  enormously  affected  by 
the  presence  of  the  wire.   We  may  think  of  the  insertion 
of  the  wire  as  resulting  in  the  transfer  of  the  e.m»f» 
f'    {E  +  -a  kB)  '  dx      to  the  ends  of  the  wire.   Of  course, 
this  formula  is  only  approximately  valid  depending  on 
the  finite  conductivity,  size,  shape,  and   velocity  of 
the  wire.   The  mathematical  idealization  of  a  perfectly 
conducting  thin  wire  is  brought  home  when  we  realize 
that  closing  the  loop  has  the  effect  of  not  allowing 
5  to  change  at  alii 
The  homo polar  generator. 

The  ideas  above  have  a  simple  application  in  the  treat- 
ment of  the  homopolar  (or  unipolar)  generator.   Typically, 
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such  a  generator  consists  of  a  perfectly  conducting  disc 
rotating  about  Its  axis  with  constant  angular  velocity 
Lj     In  the  presence  of  a  imiforin  magnetic  field  B 
parallel  to  the  axis  of  rotation.  (The  existence  of  the 
uniform  magnetic  field  extending  across  vacuum  and 
conductor  Is  consistent  with  equation  (10.7)  for  a  moving 
perfect  conductor. )  An  e»m»f«  is  produced  between  the 
axis  and  the  rim  of  the  disc. 


Figure  13»   Homo polar  generator* 

Consider  a  voltmeter  connected  by  two  perfectly 

conducting  leads  to  a  point   A  touching  the  axis  and 

a  point  B   touching  the  edge  of  the  discj  voltmeter 

and  leads  are  assumed  fixed  (Figure  13).   Let   C  be 

a  curve  which  is  carried  by  the  disc,  and  extends  from 

the  center  to  the  rim.   Since   (E  +  u  '^B).  =0  we  have 

an  e.m.f.  ,  ^yE»dx  =  -^(u'^B)«dx«   The  magnitude 

C  C 

of  the  e«m«fo  from  the  center  to  the  rim  is 
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R 

J  B  u)r  dr  =  -  (^  R^  B.   It  Is  asserted  that  this  e.m.f. 

0  '^ 

is  exactly  what  is  measured  by  the  voltmeter.   To  prove 

the  assertion  we  may  use  Faraday's  law  (10.5)  either  for 

a  moving  or  fixed  circuit.   First  let  us  take  the  path 

of  Integration  as  A  A  0  P  B  B   where  OP  designates 

the  moving  arc  C  and  PB  the  circular  arc  of  varying 

length  along  the  perimeter  of  the  disc.   The  value  of 

dj/dt  linked  by  the  circuit  made  by  this  path  completed 

via  the  voltmeter  is   B  times  the  rate  of  change  of  the 

1     2 
covered  disc  area,  ^  O  R  .   On  the  other  hand 

J  (E  +  u  < B )  •  dx  vanishes  on  each  segment  of  the  open 

path,  so  that  we  obtain  across  the  voltmeter, 

B' 

(10.8)       / ,  E  •  dx  =  ^  GJ  R^  B   . 
A  '^ 

Alternatively,  we  may  take  the  fixed  path  A  A  0  B  B 

where   OB   is  a  fixed  segment  not  moving  with  the  disc. 

Now  we  have   dj/dt  =  0.   Since   (E+u<B),  =0   on  the 

surface  of  the  disc,  we  have   E.  =  -  (ux.B),  =  WrB 

adjacent  to  the  disc  and  the  value  of 
B  B 

/  (K  +  u  X  B)  .  dx  =  /  E  •  dx  is  ^  CJ  R^  B.   Since  the 
0  0  "^ 

rest  of  the  circuit  makes  no  additional  contribution 

to  /e«  dx,   the  result  (10.8)  follows  again. 

The  electric  field  B  in  the  entire  space  is  the 

harmonic  vector  which  takes  the  given  boundary  values 

E^  =  6^;  r  B   on  the  disc. 
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The  situation  is  more  complex  when  the  magnetic 

field  B  is  nonimiform.   The  normal  component  B   is 

constant  in  time  at  a  given  moving  point  of  the  boundary 

(from  equation  (10,5) ).   Prom  this  bovmdar^  condition 

and  whatever  other  constraints  may  be  imposed  on  B, 

we  could  determine  the  magnetic  field  throughout  space 

as  a  function  of  time.   The  e.m.f.J^E'dx  =  -/(uxB)  "dx 

C         "C 
depends  only  on  the  component  B   and  hence  is  constant 

in  time  for  a  curve  C  carried  by  the  disc.   In  general, 

J  E  •  dx      depends  on  the  path  and  not  on  the  endpoints  of 
C 

C   alone.   Furthermore,  since   B   is  a  function  of  time 

the  flux  linked  by  a  given  circuit  and,  consequently, 

the  e.m.f.  registered  at  the  voltmeter  are  not  easily 

computed  in  general.   We  may,  however,  take  the  case  in 

which  the  wire   AA   is  very  close  to  the  radius   OB 

so  that  the  fluxed  linked  in  the  circuit  A  A  0  B  B 

is  negligible.   Exactly  as  In  the  prior  instance  it 

follows  that  the  voltmeter  registers  an  e.m.f.  of 

magnitude 

t 

r  E.  dx  =  oj  J  r  B^(r,t)  dr 


where   B  (r,t)   gives  the  values  of  B    along  the 
radius   OB   at  the  instant   t. 

Perhaps  it  should  be  added  that  although  the 
magnetic  field  is  "carried"  with  the  conductor  in  the 
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sense  that  the  value  of  B   follows  the  motion,  it  is 

n 

improper  to  think  of  individual  magnetic  lines  following 
the  conductor.   For  suppose  we  identify  a  given  line  by 
the  point  at  which  it  intersects  the  conductor.   It  will 
then  follow  that  as  the  given  point  on  the  disc  passes 
under  the  wire  AA  ,   the  magnetic  line  "snaps"  from 
one  side  of  the  wire  to  the  other  (see  Figure  ll|.).   In 


i  i  ^ 


nil  ii    .1   ti 


y  xxy  X  XXX  XV  y  X  yy  yx 


Figure  lii..   Magnetic  lines  near  conductors. 

other  terms,  if  a  magnetic  line  intersects  two  conductors 
it  is  clear  that  it  cannot  generally  be  assumed  to  be 
carried  by  both. 

The  argument  given  above  holds  exactly  for  a  permanent 
magnet  rotating  in  its  own  field  as  well  as  for  a  con- 
ducting disc.   It  Is  worth  remarking  also  that  perfect 
conductivity  is  not  required  for  these  arguments  provided 
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that  no  currents  flowj  we  would  still  have  the  condition 
E+u>^B=0  in  a  conductor  of  finite  conductivity. 
The  bar  magnet  is  a  case  in  point}  no  currents  will  flow 
provided  that  the  field  of  the  magnet  is  axially  symmetric , 

A   very  similar  situation  holds  if  we  have  a  rotating 
dielectric  instead  of  a  conductor.   Let  the  dielectric 
constant  be  denoted  by  K  *=  K/K  ^     where   K^  is  the 
permittivity  of  free  space*   Instead  of  Ohm's  law,  we 
have  D  =  KE   in  a  frame  of  observation  moving  with 
the   dielectric.   However,  since   D   does  not  convention- 
ally include  polarization  charge,  the  "app^iied"  field 
observed  in  the  moving  frame  is   D  =  K  u<B,   This 
leads  to  a  polarization  charge  sufficient  to  reduce  the 
electric  field  in  the  dielectric  to   E^  =  D/k  =  u k B/K. 
Returning  to  the  stationary  frame  (in  which  the  same 
polarization  charge  is  found,  non-relativistically )  we 
observe   E  =E  -u^B  =  (??  -  1)  (uxB).   For  large 
values  of  the  dielectric  constant  K,   the  result  is 
almost  the  same  as  for  a  conductor.   Of  course,  this 
must  be  measured  by  a  voltmeter  which  draws  no  current. 
The  cylindrical  pinch. 

Let  us  consider  the  problem  of  a  cylindrically 
confined  plasma  (i.e.,  a  perfectly  conducting  gas)  in 
a  cylindrical  geometry  where  the  ignorable  variable  is 
the  axial  coordinate   z.   The  plasma  cylinder  is 
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surroiinded  by  a  concentric  perfectly  conducting  shell 
(Figure  l5)»   Let  the  variable  radius  of  the  plasma  be 


Figure  l5*   Cylindrically  confined  plasma. 

denoted  by  R,   that  of  the  fixed  outer  conductor  by  a. 
Our  problem  is  to  determine   R   as  a  function  of  time 
\mder  suitably  posed  conditions.   Qualitatively,  the 
situation  is  described  as  follows.   The  vacuum  magnetic 
field  B   surro;inding  the  plasma  is  completely  determined 
by  the  external  constraints  (e.g.,  applied  e.m.f. 's 

A 

[  P'jE]      maintained  by   a  condenser  discharge)   and   the 
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plasma  radius  R,   The  magnetic  field  provides  a  driving 
pressure  B  /S^,   at  the  surface  of  the  plasma.   It  is 
intuitively  clear  that  knowledge  of  the  driving  pressure 
as  a  function  of  R  and  t   allows  one  to  compute  the 
motion  of  the  plasma  from  a  given  set  of  fluid  equations e 
A  qualitative  description  can  also  be  given  in  another 
way.   The  forces  acting  can  be  computed  as  above  using 
the  magnetic  pressure  (to  which  the  Maxwell  stress  tensor 
is  reduced  in  this  example)  or  by  using  the  formula  for 
the  force  between  current  elements.   The  "pinch"  occurs 
because  the  current  elements  in  the  plasma  c^ttract  one 
another.  However,  one  must  take  care  to  include  the 
forces  due  to  the  currents  of  the  return  circuit  in 
the  concentric  shell  and  also  balance  these  forces 
against  the  fluid  pressure  of  the  plasma  itself. 
Quantitatively,  the  first  picture  is  much  easier  to 
apply. 

In  the  vacuum  domain  D  between  the  two  cylinders 
it  is  possible  to  choose  as  a  basis  for  \J^(     ^nj  line 
P,   parallel  to  the  axis  and  any  closed  curve   P^ 
circling  the  plasma  cylinder.   For  the  open  siirface 
group  we  may  form  a  basis  of  the  surface  XIi»   "^'^^ 
annular  cross-section  of  D,   and  XIp*   ^^®  infinite 
plane  strip  joining  the  bounding  surfaces   S   and   S 
in  an  axial  section  (Figure  l6)o 
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Figure  16 «     Horaolegy  of  the   domain  between  tw®  cylinders. 

The   relevant   periods  will  be   fluxes     Jt    ~   [T"-i  3B] » 
1.2  ^   ^'^2^^^      orm.m.f.'s      [P^JB],      [P^iB].      In 
addition,    for  a  well-posed  problem   of   the   type   of 
Theorem  12   it   is   necessary   to   prescribe   Oom.f.'s 
[P^jE+u<B],      [p^iE+n<B]     where       P^      and      P^ 
are   the  boundaries   of     ^Zn       ^^*^     XIp»    respectively. 
To   avoid  infinities   in  taking  integrals  on     ^Zp      ^^^ 
P,      these  will  be   given  for   a  finite   length     "K     along 
the   axis.      Mathematically  this   corresponds   to    taking  as 
the   domain     D     the  topological   torus   defined  by  identify- 
ing corresponding  points   of  two   transverse   sections  at  a 
length     %     apart. 
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We  shall  consider  a  problem  of  the  type  of 
Theorem  12  in  the  domain  D.   On  the  boundary 
S  =  S   +  S   we  put  B^  =  0.   A  set  of  /^-periods  is 


n 

A 


given  initially  and  the  e.m.f.'s   [  Pv»E  +  u<B]   are 
prescribed  for  all  time.   Under  these  conditions,  the 
magnetic  field  between  the  plasma  and  the  outer  shell 
is  \miquely  determined.   Physically,  this  corresponds 
to  the  specification  of  the  externally  imposed  e.m.f.'s, 
for  example,  such  as  might  be  provided  by  a  condenser  dis- 
charge and  coils.   No  contribution  is  made  to  the  periods 
[  Pt^^JE  +  u><B]   by  integrals  on  the  plasma  boundary  S" 
since   (E  +  u^B),  =0  there;   the  contributions  to 
these  periods  on  the  fixed  conductor  S   will  differ 
from  zero  only  if  it  Is  slit  and  a  potential  difference 
applied  across  the  gap. 

Proceeding  from  the  existence  theorem,  the  detailed 
relations  between  the  field  and  the  given  periods  are 
easily  described.   The  time  derivative  of  the  flux  J-, 
through  y~,   is  given  by  (10.5)  as 

2-1 
=  i/u  (E+uxB)»dx-y,,  (E+u<B).dx 

where  the  boundary   P^  =  Po  -  V^p     °^     ^,±     ^^   composed 
of  the  circles   Tp   on  the  outer  conductor   S   and   Hp 
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I»  It 

on  the   plasma  boundary     S   •      On     S        we   have 

(E  +  uxB).    =  0,      and   since     u  =  0     on     S       we   obtain 
(10.9)  ii  =  -  Zi    E^  ds  =  -  Vq 

where     V^     denotes   the   e,in,f«    applied  to   the   circuit 
Pjj      aroiind   the   fixed  outer   shell   and   is    assumed  to  be 
a  known  function  of  time.      If     J?     is    the   initial  value 

of     5-      we  obtain     J,      as  a  function  of  time  by 

t 
(10.10a)  |^(t)   =1°   -JVg    {-^   dr    . 

0 
In  the  same  way,  we  may  calculate  the  time  variation 

of  the  flux  ^2  through  a  length  \     of  the  ribbon 

surface  XIp*   This  will  depend  on  the  e.m.f. 

(t)  -  J  ,    E.  ds   applied  to  a  length  "X  of  the 

Pi 

conducting  shell,  namely, 

t 

(10.10b)     i„(t)  =1°  -/v'  (t-)  dr  , 

where  \^      Is  the  initial  value  of  Jp. 

In  this  simple  geometry  the  harmonic  vector  field 
B  can  be  computed  explicitly  in  terms  of  either  fluxes 
or  currents.   The  axial  component  B   is  constant  in 


I 

V 
z 


space  and  Is  given  by 


(10.11a)     B   = 


^   Ti(a2-R2) 
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Denoting  the  azimuthal  current   about   a  length     \     of 
the   outer  conductor  by      I^     we  have   also 

(10.11b)  Bjj   =  -^     , 

whence 

(10.12)  i^  =  Ii(a2-R2)B2   =  Bt   (a2-R2)   i^     ., 

The  component  B-—  B„   =  B^  of  the  magnetic  field 

is  harmonic  in  the  plane  annular  surface  J^^      and  is 

defined  as  a  function  of  radius  r  in  terms  of  the  flux 
ig  by 

(10.13a)  Bg  =  ^— ^        , 

\r  log   (g) 

ti 

or,    in   terms   of   the   axial   current      I        in  the   plasma,   by 

z 
II 

Combining  these  results  we   obtain  the  relation 

(lo.ii^)  ia  =11  ^^°s|)  i"    . 

The  inductance  matrix  obtained  from  (10.12),  (10.14) 
is  diagonal  J 

A(a2-R2)      0 

\    °      2i  ^°S  p 

From  the  relations  (10.9),  (10.10)  we  obtain 
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(10.15)  v;  =  |^(f  (log  |)i;],  v;  =  |^t^(a2.R2)i;]  . 

If,  for  example,  both  the  current   I   and  e.in.f#   V 
are  experimentally  measured,  the  radius  of  the  plasma 
R  can  be  determined  explicitly  as  the  solution  of  this 
first  order  differential  equation.  By  employing  one 
additional  relation  which  we  shall  get  from  the  fluid 
dynamics,  we  shall  be  able  to  express  everything  in 
terms  of  a  single  type  of  quantity,  e.g.,  the  applied 
e.m.f.'s  as  known  functions  of  time. 

Thus  far  we  have  used  no  fact  about  the  plasma  except 
perfect  conductivity.   In  a  subsequent  note  we  shall  show 
that  the  quantity  p  =  p  +  3-  ,   where   p  is  the  fluid 
pressure,  is  continuous  across  the  interface   S  •   In 
particular,  we  have 

where   B   is  the  magnetic  field  at  the  plasma  side  of 

11 

the  boundary   and      B        is   the   value    at   the  boundary  of 

the    already  described  field   in     D.      Specifically,    from 
(10.11)    and    (10.13), 


(10.17) 


2TiX(a^-R^)        i^.Tt7LR^[log  a   -  log  R] 
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If  we  can  compute   p   In  terms  of  R   and  t  from  the 
fluid  side,  we  have  the  desired  relation  enabling  us  to 
solve  for  all  quantities  in  terms  of  the  applied  e.m»fo's, 
assumed  given,  or  any  other  set  of  given  data  (e.g., 
currents  )• 

Since  we  do  not  wish  here  to  consider  fluid  magnetic 
problems  in  detail,  we  proceed  by  making  simplifying 
assumptions  which  allow  us  to  obtain  ^  as  a  function 
of  R.  Vfhatever  the  nature  of  the  fluid,  it  is  clear 
in  general  that  the  complete  solution  to  the  problem 
will  be  determined  by  (10.17)  if  p   is  known  as  a 
function  of  R   and  t.   Let  us  assume  that  the  motion 
is  so  slow  that  the  fluid  density  p   is  uniform  across 
the  plasma.   Conservation  of  mass  yields 

p  =  po  ^y^^ 

where   p^,  R^   denote  initial  values  for  the  plasma. 
If  it  is  further  assumed  that  the  adiabatic  law  holds 

P  =  Pq^p/Pq^^ 

then  we  have  a  relation  between  p  and  Rj 

R 
(10.18)  p  =  Pq  {-^) 


.^O.^Y 


We  have  only  to  determine  the  plasma  boundary  field 

All 

B    to  have  a  relation  defining  R   in  terms  of  the 
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period  data  In  the   plasma  and   vacuum  domains. 

For   simplicity  we  suppose   that   the   plasma  field 
B     Is  uniform      ;    that   Is,    that  the   plasma  currents   lie 
on  the    surface      S   .      Let     £~     be   the   constant   flijix 
through  the   plasma  cross-seotlonj   hence, 

(10.19)  B  =  B^   =  — ^     . 

Putting  the  given  Information  into  (10.16)  we  obtain 

Rq  2y     i? 

(10.20)  p  =  Pq  (■#)   +72^^  • 

2ic  %lR^ 

In  a  high  temperature  plasma  y  ~  2,     and  we  obtain 


especially  simply 

(10.20a)  p  =  p^j  (^) 


^  .^0> 


where      p„      Is  the    initial   value   of     p.      Putting  the 

value   of      p      into    (10.1?)   we   obtain   a  relation  which 

defines      R     ultimately   in  terms   of    the    e.m.f.'s   in     D. 

If     B        vanishes   in   the   domain     D     between   the   plasma 
z 

and   the    outer   conductor   the   value  of     R      is   related 
directly   to    the   primary  plasma  current,   namely, 

R     2y  i^        _  i^ulf 

27.  Since    the   density    is   assumed  uniform,    this   will  be 
true    if     B     is  uniform  at    any  given   time. 

28.  The   gas   compresses    two-dlmensionally   if   the  mean 
collision   time    is   large. 
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In  case     y  ~  2,      this  equation  yields 
R  =  2%nl  /2^/l] 


II 

z 


In  any  event,    the  relation  between  the  plasma  radius 
R     and  the   currents  is   an  algebraic  one. 
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11 •      Relation  between  the   Mazwell   and  pre-Maxwell   Systems. 

It  has    already  been  remarked   that    the   type   of   data 
which   is   appropriate   to    a  well -posed  problem  for   the 
Maxwell   equations   is   entirely  different  from  that  required 
for  the   pre-Maxwell   system   (of.    Section  1).      Consequently, 
it   is  not   immediately  evident  under  what   circxmstances 
a  solution  of   the  pre-Maxwell  system  can  approximate  a 
solution  of  the   Maxwell   systemo     We    shall   show   in  a 
large  class   of  electromagnetic  problems  which  are  well- 
posed  for   the  Maxwell   system,   how  related  well-posed 

r 

problems   for   the   pre-Maxwell   system  may  be   obtained  by 
omitting  certain  of  the    specified  data  which   single   out 
a  Maxwell   solution.      In  the    sense   in  which  the   related 
pre-Maxwell   solution  is   a   satisfactory   approximation  to 
the  Maxwell    solution,    the    set  of  data  which  was   omitted 
from   the   Maxwell   problem  must   be   practically   irrelevant. 

For  convenience   we   denote    the    solution  of   a  Maxwell 
problem  by     Ej^^,   B^^,      the    solution  of    the   related  pre- 
Maxwell   problem  by     Ep,    Bp.      The    primary   question   is, 
in  what   sense    are      E      =  Ej^  -  Ep      and     B      =  Bj.  -  Bp 
small.      For    simplicity  we   take   a  fixed  bounded  domain 
D     and   assume      q=0,J=0. 

One   set   of  prescribed  data  for   a  well-posed 
Maxwell  problem  consists   of   initial    values   of     E     and 
B      subject   to      div  E  =  0      and     div  B   =  0,      and 
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tangential  boundary  values     E^^     prescribed  as   a  function 
of  time.      The   related  pre-Maxwell  problem  is   defined  as 
follows.      The   same  boundary  values     E^      are  prescribed 
for     Ep     as  for     Ej^.      It  is   easily  seen   that  the   periods 
[SjEjy,]      are   constant   in  time    (this   is   essentially  charge 
conservation)   and  the  values  of   these  periods  are 
computed  from  the   initial   values  of     E^     and  assigned 
to     Bp     for   all  time.      Prom  the   initial  values     B(x) 
we   take   only  the   boundary  values     B        and  the  periods 
[T^jB]      (which  are  well-defined  for     Bj,      as  well   as  for 
Bp     since      div  B  =  0      everywhere).      This   completes   a 
well-posed  pre-Maxwell  problem  according  to   Theorem  10. 

Another  well-posed  Maxwell   problem  requires   the 
same   initial   data  as   above,    together  with  normal 
boundary  values     E       and     B        (subject   to      [S^jB]   =  0, 
[S*}E]    =0)      and   periods      [  PjB] ,      [  PjE]      all  given 
as  fimctions   of  time.      The   closed  curve  periods   are 
defined  in   the   extended   sense    (cf.    Section  ii.)   and  must 
be    given  on   selected  curves      r'j,n'j«      The    given 
boundary  values     B        determine      E^.      to  within  an 
arbitrary  local   surface   gradient    and   the  periods 
[PjB],    [  PjE]      fix   the   periods   of   this    surface   (gradient. 
The   related  pre-Maxwell   problem  is   defined  by  prescribing 
the    same   boundary  values   and  periods      [  PjE]      and      t  PjE] 
as  in  the  Maxwell  problem,    together  with  initial  values 
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of  the  periods      [^ZJ^]      computed  from  the    initial 
values  of     B^^     (of.    Theorems  9,    11). 

We  now  state   the   following  conjectiireo      If  we 
let   the   Maxwell  problem     M     be  fixed  except    tiiat   the 
speed  of  light,      c,      is   allowed   to  become   infinite, 
the   sequence   of    solutions     RyrCc),   B^(c)      approaches 
weakly  the  limit     Ep,   Bp     of  the    associated  problemj 
(c      does  not    appear   in   the  pre-Maxwell  equations). 
More   precisely,    the   statement   is 

J       B     dt   >   0    ,  \      B     dt   >  0 

at  any  fixed  point  x  and  for  any  fixed  time  interval 

(t^.t^). 

Qualitatively,    the   above   conjecture   is   that   a 
solution  of  Maxwell's   equations   differs   from   an 
appropriately  chosen   solution  of   the   pre-Maxwell 
equations  by   a  rapidly  oscillating  function  with  mean 
zero  but   not  with  zero  mean  energy.      In  other  words, 
the    solution  of  Msixwell's   equation   separates 
unambiguously  into    a  "slowly  varying  part"   and  into 
electromagnetic   waves  which  propagate  with  the    speed 
of  light.      In  this    context,    it   is    illiominating  to  note 
that   the   solution     Bp     is    at   each   instant   the   projection 
of     B-,     onto    the    surbspace   of  harmonic   vectors    (it  is 
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easily  seen  that  the  decomposition  B^  =  B  +  Bp  is 
the  representation  of  K.     as  the  sum  of  projections 
into  the  classes  of  Table  IV  where   B   and  B 
correspond  to  the  second  and  third  columns  of  the  table 
respectively)o 

A  proof  of  the  conjecture  could  probably  be 
obtained  in  the  special  case  where  the  pre -Maxwell 
solution  is  time  Independent,  namely  for  B   independent 
of  time  (or  E.  =  0)»   By  separating  time  as  a  factor 
e       we  obtain  a  solution  of  the  Maxwell  problem  as 
an  eigenfunctlon  expansion*   The  first  terrr.»   Co  =  0, 
Is  exactly  the  pre -Maxwell  solution;  the  other  terms 
oscillate  more  and  more  rapidly  as  c  — ->  m      and 
in  the  limit  will  consequently  yield  mean  value  zero 
over  any  finite  time  Interval. 

For  an  exterior  domain,  the  conjecture  can  be 
strengthened;  the  convergence  will  be  pointwise  rather 
than  weak  since  the  initial  electromagnetic  wave  energy 
will  radiate  out  to  infinity « 

If  the  boxindary  of  the  domain  D  is  permitted  to 
move,  it  is  likely  that  the  electromagnetic  wave  component 
e' ,  b'   of  Ej,,  B^     will  vary  "adiabatically'*  with  the 
motion  of  the  walls;  loe*  the  variation  in  amplitude  of 
each  mode  can  be  computed  from  the  work  done  by  the  moving 
wall  against  the  radiation  pressure  of  the  given  wave 
component. 
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Appendl3c.   Note  on  generalization  of  the  Lagrange 
multiplier  i*ule» 

In  Section  8,  we  have  seen  how  the  transfer  of  an 
admissibility  condition  to  the  variational  side  (with 
replacement  by  the  dual  condition  according  to  Table  V) 
led  to  a  change  in  the  form  of  the  variational  fxinction 
which  may  be  described  by  the  Lagrange  multiplier  rule. 
In  general,  the  change  is  not  so  describable  and  will 
require  a  generalization  of  the  rule. 

The  simplest  cases  occur  when  the  differential 
equations  and  boundary  conditions  are  homogeneous, 
curl  X  =  0,  div  Y  =  0  and  either  X^  =  0  or  Y^^  =  0 
on  S  •  In  these  problems,   (X,Y)  may  be  written  In 
terms  of  periods  alone,  (6,5)f  (6,6).  Take  the 
admissible  category  as  U  =  Y  and  the  variational 
category  as  V  =  X.   If  periods  of  X  are  to  be 
specified,  the  term   (Y,Xq)  in  the  variational  function 
has  the  exact  form  of  the  Lagrange  multiplier  rule.   If 
periods  of  Y  are  specified,  Xq  can  be  taken  identi- 
cally zero  and  the  term   (Y,Xq)   dropped.   The  situation 
is  entirely  analogous  if  we  take  X  to  be  admissible 
and  Y  variational.   In  the  more  general  case  of 
Section  9,  we  can  still  express   (X,Y)   in  terms  of 
periods,  (9.3),  if  the  differential  relations  and 
boundary  conditions  are  homogeneous.   Those  terms  in 
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this  bilinear  sxitn  in  which  periods  of  the  variational 
class  are  given  take  the  form  of  the  Lagrange  multiplier 
rule,  while  the  tenns  which  involve  prescribed  periods 
of  the  admissible  class  are  constant  and  may  be  dropped 
from  the  variational  ftinction. 

In  the  inhomogeneous  cases,  the  more  complicated 
representations  of  the  inner  products  (6<,l5)>  (6el8) 
indicate  the  necessity  for  introducing  appropriate 
scalar  or  vector  potentials  in  order  to  express  (TJ,V) 
explicitly  in  tenns  of  the  data  which  define  the  vari- 
ational class.   If  both  5"  and  ^  differ  x'rom  zero, 
a  form  only  vaguely  reminiscent  of  the  Lagrange  multi- 
plier mile  can  be  obtained©  The  functional  may  be 
expressed  in  terms  of  certain  "partial"  potentials  ^ 
and  4  which  refer  not  to  U  and  V,   which  have  no 
potentials  in  general,  but  to  components  of  their 
orthogonal  decompositions  according  to  Table  IV,  page  73» 

If  only  one  of  C  or  "^     is  zero,  two  different 
situations  may  arise.  When  the  inhomogeneous  differ- 
ential relation  is  on  the  variational  side,  then  P[U] 
has  an  explicit  expression  as  a  functional  of  a  potential 
of  XJ   (usually  not  single-valued)^   Qither  as  FL)^]   or 
P[A],  rather  than  of  U  itself,  as  in  the  variational 
problem  of  (8«16),  (8.17),  for  example.   If,  on  the  other 
hand,  the  inhomogeneous  relation  is  in  the  admissible 
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class,  then  U  itself  is  the  argument  of  P,  but  the 
given  inhomogeneous  boundary  or  period  data  in  the 
variational  class  appear  in  the  ftmctional  as  multipliers 
which  are  defined  in  terms  of  the  potentials  of  a 
properly  chosen  variational  vector  V,  rather  than  in 
terms  of  V  itself. 
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